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1.  INTRODUCTION 


The  Flight  Dynamics  Directorate  at  Wright  Laboratory  (WL/FIBGC)  hzis  long  had 
an  interest  in  the  dynamic  and  control  aspects  of  large  space  structures  [1,  2,  3]. 
In  1985,  the  Large  Space  Structures  Technology  Program  (LSSTP)  was  initiated. 
One  of  the  goals  of  this  program  has  been  to  develop  a  laboratory  in  which  structures 
dynamically  similar  to  future  LSS  can  be  studied.  The  12-meter  truss  described  here  is 
part  of  this  effort,  involving  both  active  and  passive  methods  of  vibration  control.  WL 
has  developed  the  design  and  model  of  the  structure  and  the  :  ■‘^up  of  the  experiment. 
The  Ohio  State  University,  Control  Research  Laboratory  researchers  have  designed 
the  active  controllers,  implemented  the  softwzire  for  control,  and  performed  the  closed- 
loop  testing. 


1.1  Overview 

The  present  report  summarizes  the  total  effort  in  this  project,  including  the  theory, 
design,  simulation  studies,  and  implementation  results  of  active  vibration  control  of 
a  12  meter  cantilevered  truss  structure. 


1.2  Large  Flexible  Structures 

There  are  typically  three  areas  which  compose  Large  Space  Structure  (LSS)  control: 

•  Vibration  damping 

•  Pointing  (slewing) 

•  Static  shape  deformation  adjustment 

Each  of  the  above  requires  effort  in  modeling,  control  design,  actuation  and  sensor 
development,  and  resolution  of  practical  application  details. 

The  present  effort  addresses  the  first  area-that  of  vibration  damping. 

1.2.1  Background  and  Related  Work;  Theoretical 

Many  approaches  have  been  advocated  from  the  theoretical  viewpoint  for  vibration 
suppression  in  flexible  structures.  Most  LSS  control  schemes  for  vibration  damping 
share  in  the  goals  for  robustness  and  reduced-order  models.  During  the  last  ten  years. 
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a  number  of  approaches  have  been  suggested  and  have  been  tested  to  some  extent. 
Among  the  early  techniques  one  can  list: 

•  Basic  colocated  velocity  feedback 

•  HAC/LAC 

•  Positivity  based  feedback 

•  Independent  Modal  Control 

More  recently,  other  (mostly  state-space  based)  techniques  have  also  been  analyzed. 
The  following  should  be  highlighted: 

•  MEOP 

•  LQG/LTR 

•  Variable  Structure  (sliding-mode  control) 

•  Various  Adaptive  Control  Approaches 

•  //oo-based  techniques 

•  Decentralized  versions  of  the  above,  developed  using  different  techniques 

The  approaches  which  have  been  truly  developed  to  some  extent  for  real  systems  are 
a  small  subset  of  the  above.  In  the  present  research  we  only  had  an  opportunity  to 
consider  some  of  these  avenues. 

The  first  issue  one  has  to  face  in  designing  controllers  for  flexible  structures  is  that 
of  modeling  and  model  reduction.  In  fact,  in  a  number  of  approaches,  modeling  and 
model  reduction  are  coupled  into  the  controller  design. 

In  general,  a  Finite  Element  Model  (FEM)  is  developed,  tested,  and  provided  to  the 
control  designer.  At  this  point,  one  can  simply  truncate  the  dynamics  associated  with 
higher  frequencies.  For  beam-like  structures,  this  may  be  a  valid  approach  because 
various  criteria  based  on  control  authority,  energy  content,  etc.,  indicate  retention  of 
the  mode  ordering  sequence.  For  more  complex  structures,  however,  one  will  have 
to  utilize  an  approach  based  on  balanced  reduction  or  some  “cost”  comparison- based 
scheme.  When  considering  decentralized  control  we  used  an  approach  that  our  Group 
has  utilized  previously  at  JPL,  decentralized  balanced  realization. 

Two  other  approaches  were  also  evaluated:  Maximum  Entropy/Optimal  Projection 
(MEOP)  and  an  //oo-based  method  using  the  GAP  Metric.  In  both  of  these  cases. 
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model  reduction  and  modeling  errors  are  part  of  the  controller  design.  Modeli^.g 
errors  and  effects  of  “discarded”  modes  are  also  considered  to  an  extent  by  frequency 
shaping,  which  was  used  with  decentralized  control. 

Although  many  feedback  controllers  can  be  designed  for  decentralized  systems,  so 
far  the  most  popular  has  been  beised  on  Linear  Quadratic  cost  criteria.  Here  we 
also  have  considered  decentralized  Linear  Quadratic  control  and  utilized  a  locally 
developed  software  package  named  DOLORES.  One  aspect  considered  in  this  context 
was  frequency  weighting. 

Early  versions  of  frequency  weighting  may  be  attributed  to  the  High  Authority  Con¬ 
troller/  Low  Authority  Controller  approach  introduced  by  Aubrun  [4]  and  covered  in 
some  detail  in  the  ACOSS  Three  (Phase  I)  Final  Report.  As  analyzed  there,  Low- 
Authority  Controls  provide  limited  (broadband)  modal  damping  and  may  be  synthe¬ 
sized  using  perturbation  techniques.  It  has  been  shown  that  for  systems  restricted  to 
colocated  rate  sensing  and  decentralized  feedback  loops,  stability  is  guaranteed  over 
the  system  bandwidth.  High-Authority  Control  provides  arbitrary  damping  ratios 
in  low  frequency  modes  by  exploiting  knowledge  of  system  dynamics  and  is  gener¬ 
ally  synthesized  as  a  multi-input,  multi-output  (centralized)  control.  The  spillover 
phenomena  in  the  HAC  is  handled  by  utilizing  frequency  shaping  techniques  such  as 
presented  in  [5].  We  have  reported  preliminary  results  in  designing  the  High  Author¬ 
ity  Controller  also  assuming  decentralized  implementation  in  [6]  where  a  decentralized 
quadratic  regulator  with  frequency  weighting  was  analyzed.  As  covered  by  the  soft¬ 
ware  package  DOLORES,  various  options  can  be  analyzed  and  one  can  be  selected 
for  implementation.  Another  of  our  decentralized  control  efforts  is  beised  on  the  over¬ 
lapping  decomposition  approach-a  very  promising  idea  which  has  only  recently  been 
developed.  Overlapping  decompositions  were  first  utilized  by  Ikeda,  Siljak,  and  co¬ 
workers  (7,  8,  9,  10].  Ozgviner  and  his  students  [11]  and  Young  [12,  13]  have  extended 
the  approach  further  for  decentralized  controller  design. 

In  general,  decentralized  control  implies  that  each  control  channel  is  assigned  the  task 
of  regulation  of  a  portion  of  the  state  space.  In  the  context  of  flexib-e  structures,  each 
controller  used  in  active  vibration  damping  is  eissociated  with  a  subset  of  vibrational 
modes.  This  type  of  decomposition  can  be  performed  in  an  overlapping  manner  .so 
that  some  portion  of  the  state  space  is  assigned  to  more  than  one  channel  [14],  or 
similarly,  one  set  of  modes  is  assigned  to  more  than  one  controller. 

Once  the  local  models  are  obtained  through  the  expansion  of  overlapping  decomposi¬ 
tions,  feedback  controllers  are  designed  separately.  This  leads  to  uncoupled  feedback 
controllers  in  the  expanded  state-space,  which  are  then  “contracted”  to  the  original 
state-space  for  implementation. 

Further  details  of  the  above  mentioned  approaches  may  be  found  in  the  Chapter  2  of 
this  report. 
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1.2.2  Background  and  Related  Work:  Experimental 

Experimental  studies  on  active  vibration  damping  on  laboratory  scale  structures  have 
been  continuing  for  a  number  of  years.  Recent  surveys  of  experiments  and  facilities 
are  provided  in  [15,  16].  There  are  several  experimental  structures  and  facilities  on 
which  LSS  research  is  conducted. 

At  the  Charles  Stark  Draper  Laboratory  (CSDL),  the  Observation/Control  Spillover 
Experiment  investigated  active  vibration  damping  of  a  60-inch  cantilevered  fixed-free 
aluminum  beam  [17].  Later  experiments  have  been  performed  on  the  Flexible  Satellite 
Slew  Testbed  which  consists  of  a  rigid  central  hub  with  4  equally-spaced  flexible  arms 
extending  radially  to  a  diameter  of  9  feet.  The  entire  apparatus  is  suspended  by  a 
vertical-axis  air  bearing  [18]. 

At  the  NASA  Langley  Research  Center  (LaRC)  in  Hampton,  Virginia,  several  ex¬ 
periments  have  been  conducted.  In  early  experiments,  active  vibration  damping  of  a 
12-ft-long  free-free  uniform  aluminum  beam  [19]  and  a  7-ft  by  10-ft  flexible  grid  [20] 
ha^  been  considered.  Vibration  control  experiments  have  also  been  conducted  on  the 
20-meter,  18-bay  mini-mast  testbed  in  which  the  vertical  truss  is  cantilevered  at  the 
base  [21].  The  Spacecraft  Control  Laboratory  Experiment  (SCOLE)  investigates  the 
control  of  a  rigid  platform  with  a  10-ft  beam  connected  to  the  bottom  of  the  platform. 
At  the  opposite  end  of  the  beam  is  an  offset  reflector  with  a  diameter  of  40  inches  [22]. 
More  recent  work  is  being  conducted  on  structures  in  the  Space  Structures  Research 
Laboratory  (SSRL)  at  NASA  LaRC  [23, 24].  The  Dynamic  Scale  Model  Test  (DSMT) 
Structure  is  a  500-inch  meriform  truss  oriented  horizontally  and  cantilevered  at  one 
end.  The  CSI  Evolutionary  Model  (CEM)  is  a  52-ft  truss  structure  which  is  ori¬ 
ented  horizontally  and  entirely  suspended  by  cables.  This  structure  has  a  rib  antenna 
substructure  attached  to  one  end. 

At  the  Jet  Propulsion  Laboratory  (JPL)  of  the  California  Institute  of  Technology  in 
Pasadena,  California,  the  USAF  Astronautics  Laboratory  and  NASA  jointly  sponsor 
experimental  work  in  conjunction  with  the  Large  Flexible  Structure  Test  Facility. 
Early  vibration  suppression  and  shape  control  experiments  were  performed  on  a  16- 
ft  vertical  pinned- free  flexible  beam  [25].  Also,  in  the  Piezoelectric  Active  Member 
Experiment,  piezoceramic  material  was  used  for  sensing  and  actuation  of  a  12.5-in 
horizontal,  clamped-free  aluminum  beam  [26].  Later  vibration  suppression  work  was 
performed  on  the  Flexible  Structure  Testbed  which  consists  of  a  rigid  hub  struc¬ 
ture  with  12  equally  spaced  horizontal  ribs  projecting  out  for  a  total  diameter  of  19 
feet  [27].  More  recent  experimental  structures  include  the  Modified  Astromast,  the 
Precision  Truss,  and  a  13-bay  horizontal  truss  which  can  either  be  cantilevered  or 
suspended  at  its  midspan  [28]. 

At  the  USAF  Astronautics  Laboratory  (AFAL),  the  Large  Space  Systems  Laboratory 
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has  performed  experiments  on  the  Grid  Test  Article  which  is  a  5-ft  by  5-ft  aluminum 
grid  oriented  vertically  and  cantilevered  from  the  top  [29].  Also  at  AFAL,  the  Ad¬ 
vanced  Space  Structure  Technology  Research  Experiment  (ASTREX)  Facility  has 
been  developed.  This  structure  consists  of  an  5.5-meter  diameter  primary  mirror 
support  structure  with  three  5.2-meter  graphite  epoxy  tubes  forming  a  tripod  shape 
over  the  primary  mirror.  The  tripod  structure  is  the  secondary  mirror  support  me¬ 
tering  truss.  The  entire  apparatus  is  suspended  by  a  spherical  air  bearing  with  a 
diameter  of  48  cm.  The  air  bearing  is  sits  atop  a  4.5-meter  pedestal  allowing  the 
structure  to  rotate  about  all  three  axes  [30,  31). 

At  the  NASA  Marshall  Space  Flight  Center  (MSFC)  in  Huntsville,  Alabama,  experi¬ 
ments  are  conducted  at  the  Large  Space  Structure  Ground  Test  Facility  (LSS-GTF). 
The  Active  Control  Technique  Evaluation  for  Spacecraft  (ACES)  structure  is  a  13- 
meter  astromast  with  an  offset  antenna  with  a  diameter  of  3  meters  attached  at  the 
bottom.  The  structure  is  oriented  vertically  and  is  cantilevered  at  the  top.  A  laser, 
mirrors,  and  an  optical  detector  attached  to  the  antenna  are  used  for  pointing  con¬ 
trol  experiments.  Other  types  of  experiments  performed  on  this  structure  include 
vibration  suppression,  system  identification,  and  fault  detection  and  isolation  stud¬ 
ies  [32,  33,  34]. 

Other  experimental  facilities  and  structures  exist  at  Martin  Marietta  Denver  Aero¬ 
space,  Lockheed  Missile  and  Space  Company  (LMSC),  Harris  Corporation,  TRW,  The 
Ohio  State  University  (OSU),  the  Virginia  Polytechnic  Institute  and  State  University 
(VPI&:SU),  the  Massachusetts  Institute  of  Technology  (MIT),  the  State  University  of 
New  York  at  Buffalo  (SUNY-Buffalo),  and  Texas  A&M  University. 

The  Ohio  State  Control  Research  Laboratory  has  been  involved  in  such  -b2ised  studies 
for  about  6  years,  on  structures  both  in-house  and  at  U.  S.  Government  institutions. 
The  experimental  efforts  in  validating  control  design  approaches  for  vibration  damp¬ 
ing  of  flexible  structures  have  been  very  successful.  The  techniques  considered  and 
further  developed  as  necessary,  ranged  from  simple  Linear  Quadratic  feedback  to 
variable  structure  controllers  in  a  decentralized  configuration. 

At  the  OSU  Control  Research  Laboratory,  we  have  performed  experiments  on  vi¬ 
bration  damping  on  a  free-free  beam  with  proof-mass  actuation  [35,  36],  slewing 
for  single-link  [36,  37,  38]  and  two-link  structures  [39]  and  an  experiment  involving 
multiple  mirror  orientations  on  a  slewing  flexible  link  [39]. 

In  conjunction  with  NASA  LaRC,  we  addressed  analysis  and  design  methodologies 
for  control  and  parameter  estimation  of  large  flexible  space  structures.  The  focus 
of  the  work  was  also  on  development  and  validation  of  technology,  and  subsequent 
experimentation,  that  relates  to  the  NASA  LaRC  Spacecraft  Control  Laboratory 
Experiment  (SCOLE)  [40,  41].  The  primary  objective  of  the  research  was  threefold: 
first,  to  2issess  the  performance  of  various  methods,  primarily  from  an  adaptive  control 
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viewpoint,  in  simulation  tests  and  actual  experimentation;  and  second,  to  determine 
the  effectiveness  of  various  techniques  of  parameter  estimation  for  use  in  control 
applications.  Finally,  simplified  minimum-time  slewing  for  flexible  structures  was 
also  studied. 

At  JPL,  we  investigated  various  decentralized  control  approaches  primarily  related 
to  the  overlapping  decomposition  approach.  The  main  thrust  of  the  project  was 
verification  of  the  decentralized  control  approaches  on  the  JPL/AFAL  Flexible  Struc¬ 
ture  Testbed.  Also,  some  adaptive  control  and  fault  detection  schemes  were  investi¬ 
gated  [42]. 


1.3  Description  of  the  Experimental  Facility 

The  structure  considered  is  a  12-meter  truss  constructed  from  square  aluminum  tub¬ 
ing,  with  cross  members  of  Lexan  plastic.  It  is  composed  of  four  sections-each  con¬ 
sisting  of  4  truss  bays.  These  sections  are  bolted  together  to  form  the  truss  with  a 
total  of  16  bays.  For  closed-loop  control,  the  truss  is  oriented  in  the  vertical  cantilever 
position  [43,  44). 

Active  damping  is  achieved  through  eight  momentum  exchange  proof-mass  actuators. 
These  devices  are  essentially  linear  DC  motors-each  supplied  by  its  own  current  driver 
and  capable  of  supplying  up  to  1  pound  force  from  2  Hz  to  100  Hz.  Pairs  of  actuators 
are  mounted  perpendicular  to  each  other  at  the  center  of  both  the  one-half  and 
three-quarter  stations  on  the  truss,  thus  affecting  only  the  bending  modes.  There 
are  four  actuators  at  the  tip,  two  in  each  direction,  mounted  off  center.  In  this 
way,  both  torsion  and  bending  can  be  affected.  Sensors  for  the  control  system  are 
accelerometers,  which  are  mounted  on  the  actuator  housings.  The  acceleration  signals 
are  integrated  to  provide  velocity  signals.  A  point  light  source  is  located  at  the  tip. 
An  optical  sensor  can  record  the  position  of  the  light  source  to  provide  information 
on  the  displacement  of  the  tip. 

The  truss  is  excited  via  two  alternative  sources.  The  first  is  an  additional  proof  mass 
actuator,  identical  to  the  control  actuators.  The  second  alternative  is  a  large  shaker, 
mounted  to  the  wall  at  the  tip  of  the  truss.  The  shaker  applies  a  force  up  to  10  lb 
through  a  light  connecting  arm. 

The  control  algorithms  are  implemented  digitally  using  a  custom  control  computer 
manufactured  by  Systolic  Systems,  Inc.  The  system  consists  of  the  control  computer 
and  a  Sun  3/50  graphics  workstation  that  serves  as  a  development  system.  External 
control  is  accomplished  through  a  12-channeI,  16-bit,  A/D  and  D/A  card.  The  control 
programs  are  written  in  C  programming  language,  compiled  on  the  Sun  workstation, 
and  downloaded  to  the  control  computer.  This  system  is  capable  of  implementing 
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rather  complex  control  algorithms  at  frequencies  needed  for  the  truss. 

Data  acquisition  and  analysis  can  be  performed  either  by  uploading  data  to  the  Sun 
host  or  through  an  OnoSokki  dual  channel  spectrum  analyzer.  Since  the  interface 
between  the  control  computer  and  the  workstation  is  very  slow,  the  spectrum  analyzer 
is  frequently  used  to  minimize  the  time  required  to  collect  data  for  each  experiment. 
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2.  ANALYSIS  AND  DESIGN  OF  CONTROLLERS  FOR  FLEXIBLE 

STRUCTURES 


Partial  differential  equation  (PDE)  models  and  finite  element  models  (FEMs)  are 
often  used  to  model  flexible  structures.  For  controller  design,  the  FEM  is  often  trans¬ 
formed  from  nodal  coordinates  to  modal  coordinates  by  simultaneously  diagonalizing 
the  mass  and  stiffness  matrices.  The  model,  put  into  state  space  form,  can  then 
be  reduced  using  balanced  model  reduction  techniques.  Controllers  are  subsequently 
designed  for  the  reduced-order  model. 

2.1  Modeling  Issues 

2.1.1  Partial  Differential  Equation  Model 

The  first  and  probably  most  obvious  problem  is  modeling  the  structural  dynamics. 
Large  space  structures  are  infinite  dimensional.  Using  the  notation  of  [45],  they  may 
be  described,  in  general,  by  a  system  of  partial  differential  equations  (PDE): 

m{x)u{x,t)  +  Dou{x^t)  -h  Aou{x,t)  =  F{xyt) 

where  u(x,f)  represents  the  instantaneous  displacement  of  the  structure  off  of  the 
equilibrium  position;  m(a:)  is  the  mass  density;  DqU  is  the  damping  term;  AqU  is  the 
internal  restoring  force  term;  and  F(i,f)  is  the  external  force  distribution. 

Aq  is  a  positive  semi-definite  differential  operator  on  the  infinite  dimensional  Hilbert 
space  of  square  integrable  functions  (with  the  usual  inner  product  and  associated 
norm)  defined  on  (the  interval  containing)  the  structure,  fl  :  Ho  =  The 

operator  is  usually  assumed  to  have  a  discrete  spectrum  such  that  it  satisfies  the 
eigenproblem 


Ao<Pk  = 


where  ujk  represents  the  modal  frequencies  and  <f>k  represents  the  mode  shapes. 

The  external  force  distribution  can  be  decomposed  into  control  forces  and  external 
disturbance  forces  (including  nonlinearities): 

F{x,  f )  =  Fc{x,  t)  -I-  Frf(x,  t) 

The  control  forces  are  due  to  M  actuators  and  may  be  described  by: 
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F,  =  Bof. 


where  /  represents  the  actuator  amplitudes  and  Bq  represents  the  axtuator  influence 
functions  (which,  in  general,  allow  for  both  point  and  distributed  actuation). 

Meaisurements  from  P  position  or  velocity  sensors  (including  integrated  accelerometer 
measurements)  can  be  described  as 

y  =  Cqu  4-  Eoit 

where  Co  represents  the  influence  functions  of  the  position  sensors  and  Eq  represents 
the  influence  functions  of  the  velocity  sensors  (which  also  allow  for  both  point  and 
distributed  measurement). 

State  Space  Formulation  from  the  PDE  Model 

Since  the  truss  can  basically  be  treated  as  a  cantilevered  beam,  we  shall  generate 
a  state  space  representation  of  the  system  dynamics  from  the  Euler-Bernoulli  PDE 
model.  Consider, 

(2.1) 

pA  ox*  oi^ 


where 


p  =  Density  of  Beam 
A  =  Cross  Sectional  Area  of  Beam 
E  =  Young’s  Modulus 
I  =  Area  Moment  of  Inertia  of  Beam 
w  =  Natural  Frequency  (rad/sec) 

The  solution  can  now  be  obtained  by  employing  separation  of  variables: 

u(i,0  =  r^(i)g*(t)  (2.2) 


Substitution  of  Equation  (2.2)  into  Equation  (2.1),  yields  a  time  dependent  and  a 
position  dependent  equation. 


Az(0 

df^ 


+  -  0 


(2.3) 
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(2.4) 


d*r,(x) 

dx* 


pA 

Ely 


(x)  =  0 


A  general  solution  of  the  position  dependent  equation  is  chosen  to  be 

T^{x)  =  Ax  sin  ^xX  +  Bx  cos  /3,x  +  C*  sinh  /3*i  +  Dx  cosh  (2.5) 

Applying  boundary  conditions  will  allow  us  to  solve  for  the  coefficients  and  the  vi¬ 
brational  frequencies. 

The  beam  can  now  be  forced  with  a  set  of  point  force  actuators  applied  at  specified 
locations.  By  introducing  forces,  Fait),  acting  transversely  on  the  beam  for  the  z  —  x 
bending,  at  actuator  locations  Xa  (a  =  1,2, . . .  ,Na),  and  equating  these  forces  with 
Equation  (2.1)  results  in 


El  &*u(x,t)  d^u(x,t) 

pA  dx*  dO 


-  x.)F.(t) 

P^  a=l 


where 


•=i 


Substituting  Equation  (2.7)  into  Equation  (2.6)  gives 


Ely 

pA 


dx* 


-]  <>.i( 


Na 


' d*r  CtA  °°  1 

“ ’■•iU  1  .  .(,)  + 

P^  0  =  1 


i=l 


From  Equation  (2.4), 

=  P*xrx{x) 


dx* 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


where 


Upon  substitution  of  Equation  (2.9)  into  Equation  (2.8),  one  obtains 
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For  the  orthonormal  with  respect  to 


r,i{x)r,jix)dx  =  I  j’  *  J  j  (2.1 1) 

Equation  (2.11)  is  multiplied  by  rjj(x)  and  integrated  over  f^(‘)dx  to  obtain 

|£  »'*i(a:)r.i(a:)(9,.(t)  +  u;,?9xi(<))|  = 

j5Z»‘*i(^)^(a:-®a)F„(o|<ix,  j  =  l . oo  (2.12) 

which  reduces  to 


1  Sa 

qz,{t)  +  <^zhzi{t)  =  74X1  r,i(x«)F«(<)  (2.13) 

Expanding  Equation  (2.13)  for  the  first  n-modes, 

9*1  =  ^r^i{xi)Fi(t)  +  rti{x2)F2{t)  + ’•■  +  rti{xNa)FNa{t) 

qz3  +  =  ^^*2(*l)^l(0  +  r^2{x2)F2{t)  +  • '  •  +  r,2(xyv„)Fjv„(t) 

9*„+W*n9*n  =  ^»'*n(*i)Fi(0  +  »'*n(a;2)F2(0  +  -"+r,„(xA?„)Flv.(t)  (2.14) 

and  dropping  the  z  subscript,  Equation  (2.14)  can  be  more  generally  represented  by 


L 


(2.15) 


Note  that  all  mode  shapes  are  to  be  evaluated  at  the  point  of  actuation.  The  state- 
space  model  is  therefore 


(2.16) 


n  =  i4r2  -j" 


0  1 


F 


where 


n  =  [<7  9]^ 

Ai 

Inclusion  of  Actuators  and  Sensors  into  an  Expanded  State-Space  Model 

The  sensor  dynamics  are  implicitly  included  in  the  overall  state-space  model  by  as¬ 
suming  that  the  output  of  the  colocated  sensors  is  velocity,  and  not  acceleration. 
The  sensor  scale  factors,  in  volts/m/sec,  are  the  only  values  taken  into  account,  since 
the  sensors  had  flat  frequency  responses  over  the  range  of  interest.  Only  the  output 
C,  D-matrices  need  to  be  adjusted  to  reflect  this  fact  (i.e.,  C  hcis  nonzero  elements 
corresponding  to  velocity  states,  D  is  the  null  matrix).  Therefore,  it  is  the  actuator 
dynamics  that  must  be  combined  with  the  truss  dynamics. 

Define  the  beam/truss  dynamics  to  be  the  following: 


0  /] 
-wl  0 


.Yy  =  AyXy^ByFy  (2.17a) 

.Vx  =  AxXx^BxFx  (2.17b) 

Xf  =  AjXj  +  BjFt  (2.17c) 

the  actuators  to  have  dynamics, 

Xi  =  AiXi  BiVi-Y  GyxXy  +  GTiXr;  *  =  y  direction  (2.18a) 

Xi  =  AiXi  +  BiVi  +  GziXx  +  GjiXT]  i=  x  direction  (2.18b) 

and  the  forces  acting  in  each  degree  of  freedom, 

Fy  =  CiXi  +  DiVi  (2.19a) 

Fx  =  CiXi  +  DiVi  (2.19b) 

Ft  =  CiXi  +  DiVi  (2.19c) 

where 


0 

1 

-2^w 
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^•=  X  4  *  =  1.2, 3, 4 

Q.  i  =  l,2,3,4 

Ci  =  mj,[-ki  -di]  i  =  l,2,3,4 
Di  =  rripGi  t  =  1,2,3,4 

where 

State  Vector  =  \Xy  Xx  Xj  X\  . Xs]^ 

Input  Vector  =  [Vj  . Vs]^ 

Output  Vector  =  [Y  Z  6^ 

where 

Vi  =  Voltage  input  to  power  amplifier  -  i 
2.1.2  Finite  Element  Model 

Note  that  the  description  discussed  in  the  previous  section  is  for  infinite-dimensional 
distributed  parameter  systems.  Obviously,  it  would  be  impractical  (if  not  impossible) 
to  design  controllers  for  LSS  using  an  infinite-dimensional  model.  Therefore,  a  finite¬ 
dimensional  model  is  desired.  Perhaps  the  most  popular  method  is  finite  element 
modeling  (FEM)  where  the  resulting  model  is  of  the  form 

Mz-\-  Dz^r  Kz  =  Q  (2.20) 

where  A/,  D,  and  K  are  the  mass,  damping,  and  stiffness  matrices,  respectively,  z  is 
the  nodal  coordinate  vector  (displacements  and  velocities  in  each  direction,  at  each 
node);  and  Q  is  the  generalized  force  vector.  Note  that  D  is  usually  positive  semi- 
dehnite  and  symmetric  representing  the  internal  structural  damping,  but  can  also 
contain  a  skew  symmetric  component  representing  gyroscopic  damping  due  to  any  on¬ 
board  rotors  or  constant  spin  rate  of  the  whole  LSS.  A/  is  a  positive  definite,  symmetric 
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matrix,  and  A"  is  a  positive  semi-definite,  symmetric  matrix.  These  matrices  can  be 
generated  from  any  of  several  existing  FEM  software  packages  such  as  NASTRAN, 
SPAR,  STARDYN,  etc.  Whatever  the  case,  the  structural  dynamicist  obtains  a  FEM 
of  the  form  above. 


2.2  Model  Order  Reduction 


2.2.1  Transformation  from  Nodal  Coordinates  to  Modal  Coordinates 

One  basic  approach  to  order  reduction  is  through  consideration  of  the  modal  distri¬ 
bution.  In  order  to  use  this  approach  for  the  FEM,  the  model  must  be  converted  from 
nodal  coordinates  to  modal  coordinates  by  simultaneously  diagonalizing  the  mciss  and 
stiffness  matrices.  These  matrices  may  be  diagonalized  in  a  variety  of  ways  resulting 
in  differently  balanced  relations  between  the  inputs  to  states  and  states  to  outputs 
transformations.  The  balanced  realization  obtained  here  results  in  Cpt  =  where 
Cpt  describes  the  relation  between  the  nodal  state  vector  and  a  velocity  sensor,  and 
Bpt  describes  the  relation  between  a  point  force  input  at  the  same  location  and  the 
nodal  states.  This  provides  some  symmetry  in  the  model  and  makes  the  model  more 
numerically  stable  for  simulations. 

The  state  vector  is  defined  to  be 

X  = 

Xm 

where  Xm  represents  the  modal  displacements  and  represents  the  modal  velocities. 
To  diagonalize  the  system,  the  matrix  $  is  used  to  define  the  unitary  transformation 
z  =  ^Xm  such  that 

=  I  (2.22) 


(2.21) 


and 


=  (i^  =  diag{iol  (2.23) 

where  fl  represents  the  modal  frequencies  and  n  is  the  number  of  modes  in  the  model. 

This  matrix  $  can  be  obtained  using  singular  value  decomposition  (SVD)  methods. 
First,  find  diagonal  fli  and  unitary  U\  (  (/f*  =  )  such  that 


fi,  =UlMUi. 


(2.24) 
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I,rt  U2  =  ^  and  k  =  KU1U2  such  that  K  is  symmetric  and  positive 

semidefinite.  Then,  find  diagonal  ^2  and  unitary  U3  such  that 

^2  =  UlkU3.  (2.25) 

By  defining  $  =  U1U2U3,  it  is  evident  that  (2.22)  and  (2.23)  hold.  Note  that 
is  merely  fi2-  Now  Equation  (2.20)  can  be  transformed  to  modal  coordinates  by 
substituting  z  =  and  premultiplying  by 

=  ^^BgU  (2.26) 

where  Bg  is  the  force  distribution  matrix.  Now,  Equation  (2.26)  can  be  rewritten  in 
state  space  form  using  Equation  (2.21) 

.  im  J  “  [s-n^  J  [  J  [  ^'^Bg  J  “  -  (2-27) 

The  output  matrix  for  velocities  at  the  sensor  locations  can  be  written 

im 

where  C,  is  the  sensor  distribution  matrix.  The  displacements  at  the  actuator  loca¬ 
tions  are, 

y  =  Cjx=[Cg^  0  1  "I""*  (2.29) 

Thus,  the  structure  is  represented  by  a  finite-dimensional  model  in  modal  coordinates. 
2.2.2  Balanced  Reduction 

The  modal  transformation  described  in  the  last  section  results  in  a  system  equal  in 
order  to  the  FEM  model.  This  model  has  too  many  modes  to  be  used  for  controller 
design,  and  the  higher  order  modes  are  certainly  inaccurate.  In  modal  coordinates, 
modes  may  be  neglected  simply  by  deleting  the  rows  and  columns  of  the  state  space 
representation  which  represent  these  modes.  Modes  were  deleted  according  the  their 
controllability  and/or  observability.  Since  the  actuators  and  sensors  are  colocated, 
the  controllability  and  observability  for  structural  modes  are  closely  related. 


(2.28) 
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For  the  structural  modes  (as  opposed  to  actuator  modes),  the  controllability  and 
observability  are  very  strongly  linked  to  the  modal  frequency — lower  frequency  modes 
are  most  important.  In  fact,  for  a  white  noise  disturbance  at  most  points  on  the  truss, 
the  first  few  modes  very  accurately  describe  the  motion  of  the  truss. 

The  decision  of  which  modes  to  retain  from  each  subsystem  can  be  made  by  con¬ 
sidering  the  grammians  [46].  The  basic  idea  is  to  transform  a  system  so  that  its 
controllability  and  observability  grammians  are  equal  and  diagonal  (the  system  is 
internally  balanced).  For  a  system  in  this  form,  the  singular  values  of  the  grammi¬ 
ans  are  contained  in  a  diagonal  matrix  with  the  values  in  descending  order.  With 
the  controllability  and  observability  grammians  equal,  we  can  then  make  meaningful 
statements  about  how  controllable  and  observable  a  given  mode  is,  and  delete  those 
modes  from  the  model  that  are  weakly  controllable  (and  thus  weakly  observable). 

In  general,  the  controllability  grammian  Wc  is  the  solution  of  the  algebraic  Lyapunov 
equation; 

AWc  +  A  =  0  (2.30) 

while  the  observability  grammian  Wo  satisfies  the  equation 

A^Wo  -f  Wo  A  +  C'^C  =  0  (2.31 } 


The  block  diagonal  form  of  a  system  in  modal  coordinates  (see  equation  (2.52a)  for 
example)  is  exploited  to  give  closed  form  solutions  for  Wc  and  Wg.  For  W^  the  2x2 
blocks  Wij  are  a  solution  of 

A^W,,  -f  WijAj  +  BiBj  =  0  (2.32) 

with  Ai,  Bi  in  modal  coordinate  form,  and  the  closed  form  solution  is 


/  2u)iUJj{(ju;i  +  (iUj) 
\  -Wi(u;J-a;?) 


-  u;f) 

2u;<u;j(C,Wi  -|- 


(2.33) 


where  =  bj^bj,  dij  is  a  measure  of  how  closely  correlated  modes  i  and  j  are, 

and  represents  the  modal  frequency  and  damping  of  the  ith  mode.  For  the 

observability  grammian  Wg,  Equation  (2.32)  takes  a  form  similar  to  Equation  (2.31) 
and  0ij  in  Equation  (2.33)  is  replaced  with  'fij  =  cjcj. 

For  the  case  where  p  >  m  (more  sensors  than  actuators),  it  can  be  shown  that  the 
system  is  orthogonally  symmetric  if  we  can  find  U  such  that  C  =  U B^P  and  P  is 
diag(l,-l,...,l,-l)  which  satisfies  A'^  =  PAP.  The  cross-grammian  is  now  defined 
as  Wco  which  satisfies 
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AWco  +  WcoA^  +  BU'^C  =  0 


(2.34) 


There  are  several  relations  of  Wco  that  make  it  extremely  useful  for  balancing  appli¬ 
cations: 

Wl  = 

and  Wco  =  WcP  =  PWo 

Thus  the  equation  above  for  Wij  works  for  all  three  grammians  with  appropriate 
changes  in  the  f^ij  term,  where  for  Wco,  0ij  =  7«j- 

Now  we  consider  the  balancing  problem  for  a  system  in  modal  coordinates.  The  state 
transformation  T  to  take  the  system  described  in  state  space  as  {j4,  B,  C}  to  balanced 
form  {T~^AT,T~^B,CT},  with  equal  and  diagonal  grammians 


T-^WcT-'^  =  T'^WoT  =  E  =  diag(ai) 


(2.36) 


and  with  <ti  >  <T2  ^  >  0  consists  of  right  eigenvectors  of  WcWo,  since 

I?  =  {t-^  WcT-'^)  (t'^WoT)  =  T'*  {WcWo)  T  (2.37) 

For  orthogonally  symmetric  systems,  W^^  =  WcWo  implies  that  all  balancing  infor¬ 
mation  is  in  the  eigenstructure  of  the  cross-grammian: 

T-^WcoT  =  h  (2.38) 


(2.39) 


so  A  =  diag(±<7i)  with  E  =  abs(A). 

The  above  results  can  be  used  for  model  reduction  by  choosing  the  most  dominant 
modes  by  observing  large  magnitude  changes  of  the  singular  values  on  the  diagonal. 
The  dominant  balanced  subsystem  h<is  state-space  representation 

{^.5,6}  =  |(/,0)T-Mr^  '  j,(/,0)r-'B,cr(  J  j|  (2.40) 

wliere  columns  of  T  f  ^  J  =  Ii  are  right  eigenvectors  of  Wee  corresponding  to  the 
largest  eigenvalues  of  Ai. 
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Actuator  modes  are  slightly  more  complicated.  When  the  model  is  transformed  to 
modal  coordinates,  the  actuator  modes  do  not  correspond  to  single  actuators,  but 
rather  to  a  linear  combination  of  all  actuators. 


2.3  Controller  Design  Approaches 

This  project  presented  the  opportunity  to  test  a  wide  variety  of  control  techniques. 
These  were  all  based  on  the  Linear  Quadratic  Regulator  (LQR).  Linear  quadratic 
regulator  (LQR)  design  methods  are  common  and  convenient  controller  design  tech¬ 
niques  [47,  48).  The  cost  function  minimized  in  these  techniques  is  typically 

(^x^Qx  +  Ru'j  dt  (2.41) 

or  some  variation  thereof.  Controllers  are  often  designed  using  the  pointing  require¬ 
ments.  The  inputs  are  usually  weighted  equally,  thus  R  is  diagonal. 

One  particular  method  of  interest  is  linear  optimal  output  feedback  where  only  the 
output  is  available  for  feedback  [49].  Often  in  a  LQR  design,  the  outputs  of  a  system 
are  weighted  rather  than  the  states.  The  resulting  cost  criterion  is  of  the  form 

Jo  u^Ru)  dt  (2.42) 

where  u  and  y  are  the  system  input  and  output  vectors,  respectively.  Q  and  R 
are  positive  semi-definite  and  positive  definite  matrices,  respectively,  of  appropriate 
dimensions.  If  Q  is  positive  definite  and  the  system  is  observable,  this  cost  criterion  is 
a  Lyapunov  function  and  the  resulting  controller  asymptotically  stabilizes  the  system. 

All  of  these  design  techniques  have  been  discussed  thoroughly  in  other  sources;  there¬ 
fore,  only  a  short  description  will  be  given  here. 


2.3.1  Linear  Quadratic  Gaussian  (LQG) 

In  this  approach,  it  is  assumed  that  the  system  equations  are 

X  =  Ax  -h  Bu  +  vj  (2.43) 

y  =  Cx  ->rv  (2.44) 

where  w  ~  A^(0,  W)  and  v  ~  ^(0,  V)  are  independent  white  Gaussian  processes.  The 
cost  criterion  which  is  to  be  minimized  is 
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(2.45) 


J  =  lir^ 

This  approach  uses  the  standard  full  state  feedback  with  a  Kalman  filter  to  estimate 
the  modal  states  [50].  The  Kalman  filter  is  of  the  form 

Ax -\-Bu  +  L{y-  Cx)  (2.46) 

where  x  is  the  estimated  state  vector.  The  gain  L  is  found  by  solving  the  Algebraic 
Riccati  equation, 

PA'^  +  AP  AW-  PC'^V-'  CP  =  0  (2.47) 

where  W  is  the  process  noise  covariance  matrix,  V  is  the  meaisurement  noise  covari¬ 
ance  matrix  and  L  =  PC^V~^. 

The  state  feedback  was  computed  to  minimize  the  cost  function  by  solving  a  similar 
Riccati  equation, 

SAA  A'^S  -  SBR-^B'^S  +  Q  =  0  (2.48) 

The  feedback  is  given  by  u  =  Kx  =  —R'^B^Sx. 


Linear  Quadratic  Gaussian  with  Loop  TVansfer  Recovery  (LQG/LTR) 

The  LQG/LTR  approach  [51,  52,  11]  combines  the  frequency  domain  and  state  space 
techniques  for  a  minimum  phase  system  and  is  very  similar  to  the  LQG  design  dis¬ 
cussed  above.  The  fundamental  idea  in  LTR  design  is  to  recover  a  target  feedback 
loop  (TFL)  which  satisfies  some  design  performance  specification  with  a  suitable 
asymptotic  design.  An  example  of  a  TFL  is  a  LQR  design.  If  system  states  are  not 
available  and  they  have  to  be  estimated,  some  of  the  robustness  properties  of  the 
LQR  design  will  be  lost.  LTR  provides  a  remedy  for  this.  The  theory  states  that  a 
filter  (not  a  Kalman  filter)  can  be  designed  such  that  the  overall  design — controller 
plus  filter — will  asymptotically  behave  as  the  controller  design  (LQR). 

The  primary  difference  between  the  LTR  technique  and  the  LQG  technique  is  the 
filter  design  equation 

PA'^  A  APA  BB'^  -  PC'^CP  =  0  (2.49) 

The  theory  states  that  ais  the  design  parameter  q  —*  oo,  the  system  feedback  loop 
approaches  the  TFL  asymptotically. 
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2.3.2  Decentralized  Control 


In  the  control  of  Large  Flexible  Space  Structures,  the  notion  of  decentralization  plays 
an  important  role.  This  notion  is  concerned  with  imposing  information  flow  con¬ 
straints  on  the  various  controllers  of  the  system;  controllers  of  the  system  are  only 
allowed  to  measure  certain  outputs  of  the  system  and  control  through  corresponding 
inputs.  The  decentralization  constraint  enters  into  large  scale  systems  because  it  may 
be  impractical  or  even  impossible  to  communicate  signals  from  one  controller  to  an¬ 
other.  Moreover,  decentralization  may  be  required  by  the  control  designer  to  achieve 
reliability  and  a  degree  of  redundancy,  and  it  may  impose  a  structure  to  the  control 
implementation  by  relegating  control  authority  to  separate  channels. 

A  large  scale  system  is  called  decentralized  when  either  due  to  physical  measure¬ 
ment/actuation  constraints,  or  by  design,  there  exists  a  one  to  one  correspondence 
between  sets  of  inputs  and  outputs.  If  we  assume  that  the  large  scale  system  is  linear, 
the  state  space  model  is, 

X  =  Ax  -|-  BiUi 

yi  =  CiX  -\-  DiUi  i  =  1 , . . . , 

where  the  index  t  denotes  an  input/output  channel. 

During  the  last  10  years,  a  number  of  techniques  have  been  considered  for  designing 
decentralized  controllers.  A  number  of  optimal  control  approaches  have  been  pursued 
by  Davison  [53,  54],  Ozguner  [55,  56],  and  others  [57].  These  basically  define  a 
Linear  Quadratic  (LQ)  framework  and  may  utilize  various  avenues  for  solving  for  the 
“optimal”  feedback  gains.  Servo-compensators  for  tracking  and  disturbance  rejection, 
implementation  constraints  and  frequency  weighting  may  be  included. 

Recently,  adaptive  control  techniques  have  also  been  used  for  decentralized  systems. 
Applications  of  model-reference  adaptive  control  have  been  made  to  robotic  systems 
and  flexible  structures.  Indirect  adaptive  control  has  also  been  considered  on  a  flexible 
structure. 

Decentralization  is  particularly  relevant  in  Large  Flexible  Space  Structures  for  various 
reasons; 

•  The  sensor/actuator  locations  on  the  structure  have  a  decentralized  nature 
which  makes  it  difficult  or  even  impossible  to  incorporate  a  centralized  con¬ 
trol  architecture. 

•  Information  transfer  will  be  costly  for  complex  space  systems,  so  that  the  de¬ 
signer  may  wish  to  impose  decentralized  (local)  feedback. 


(2.50) 
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•  On-line  computations  will  be  too  time  consuming  for  some  structures  so  simpler 
decentralized  strategies  will  be  necessary. 

A  model  for  the  large  flexible  structure  to  be  controlled  under  a  decentralized  infor¬ 
mation  structure  is  the  following: 

Vi  —  9  “b  Cvi  9  *  —  1 , . . . , 


where  u,-  6  1R"*‘  is  the  vector  of  control  inputs  available  at  the  channel,  y,  €  IR’^* 
is  the  vector  of  measurements  available  at  the  t*^  channel,  and  q  €  IR"  is  the  nodal 
displacement  vector.  M  and  K  are  the  inertia  and  stiffness  matrices,  respectively. 
Such  a  formulation  is  general  enough  to  include  the  cases  of  position  measurements 
only  (through  Cr,,)  or  velocity  measurements  only  (through  C„,),  as  well  as  the  case 
of  distributed  control  with  noncolocated  sensing.  Note  that  in  all  of  the  modeling 
and  control  philosophies  mentioned  above,  noncolocated  feedback  control  schemes 
may  also  be  used.  That  is,  decentralized  control  does  not  necessarily  imply  colocated 
feedback,  but  rather  the  assignment  of  one  or  more  sets  of  sensors  to  correspond  to 
one  or  more  sets  of  actuators.  The  primary  claim  with  the  formulation  (2.51)  is  that 
certain  measurements  correspond  to  certain  control  inputs  for  purposes  of  feedback 
control.  This  decentralized  feedback  mechanism  may  be  due  to  spatial  decentraliza¬ 
tion  in  the  large  structure,  or  it  may  be  due  to  an  imposed  structure  dictated  by  the 
designer,  so  that  the  channel-by-channel  feedback  is  simpler  to  implement  and  faster 
to  compute  in  real  time. 

The  well  known  transformation  q  =  where  $  is  a  matrix  of  modal  vectors  (the 
mode  shape  matrix),  is  now  utilized  as  discussed  in  Section  2.2.1  to  simultaneously 
diagonalize  M  and  K  in  (2.51)  to  give  the  state  space  formulation 


dt  c 


y. 


0  /  ■ 

C' 

0 

.C. 

C 

.  C 

N 


0 


«• 


(2.52a) 

(2.52b) 


We  note  that  in  the  formulation  (2.52a)  no  damping  is  present  in  the  structure;  for 
our  purposes,  a  term  to  account  for  damping  effects  may  be  easily  included  in  place 
of  the  lower-right  zero  of  the  first  matrix  in  (2.52a).  Typically,  damping  ratios  on 
the  order  of  0.5%  are  assumed  as  the  damping  effects  are  very  small  for  large  flexible 
structures.  Furthermore,  actuator  and  sensor  dynamics  can  also  be  appended  to  this 
formulation. 
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2.3.3  Centralized  and  Decentralized  Optimal  Output  Feedback 

Consider  the  basic  decentralized  quadratic  regulator  problem: 

X  =  Ax  BiUi  x(0)  =  xo  (2.53) 

«=i 

Vi  =  CiX-\-DiUi  (2.54) 

with  cost 

J  =  /  {x^Qx  +  2^  x^NiUi  +  Yl  ujRiUi)dt  (2.55) 

i=l  1=1 

and  the  following  feedback  structure  constraint: 

=  I<i  Vi  t  =  1 , . . . ,  I/.  (2.56) 

Remark  1;  In  the  following  analysis,  the  terms  corresponding  to  a  direct  feedthrough 
of  the  inputs  to  the  outputs  (i.e.,  Di)  are  neglected  for  notational  convenience.  If 
direct  feedthroughs  are  present,  outputs  of  the  forms  j/"'*"  =  y,  —  D,u,  can  be  utilized. 

It  can  be  shown  that  the  necessary  conditions  for  minimizing  J  given  by  (2.55)  with 
the  controller  structure  (2.56)  imply  the  solution  of  the  following  system  of  nonlinear 
algebraic  equations: 


AjP  +  PAc  +  Q  =  0 
AcL  +  LAJ  +  Xq  =  0 


(2.57) 


and 


Vfi.J  =  Bj  PLCj  +  RiKiCiLCj  +  NJlcJ  =  0  i  =  1 , . . . , »/  (2.58) 

where 

Ac  =  A-h^  BiKiCi  (2.59) 

i=l 

Q  =  Q  +  +  +  (2.60) 

i=l  i=l 

Xq  =  Xo  Xq  .  (2.61) 

A  software  package  has  been  developed  recently  [58]  to  solve  this  problem  and  has 
been  used  in  all  the  numerical  examples  that  follow.  Many  approaches  have  been 
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introduced  in  the  literature  to  solve  the  above  system  of  equations.  The  approach 
pursued  in  our  software  package  is  a  gradient  method  which  gives  at  iteration  j 

=  a:/ -e  (2.62) 

where  f  is  a  step  size  chosen  by  appealing  to  the  Fibonacci  linear  search  technique. 
The  solution  methodology  is  as  follows; 

i)  Choose  an  initial  stabilizing  gain  and  a  tolerance  6  on  the  norm  of  the  gradient 

Vk,J. 

ii)  Solve  the  Lyapunov  equations  given  by  (2.57). 

iii)  Find  the  new  gain  using  (2.62)  and  the  Fibonacci  linear  search. 

iv)  Go  back  to  step  (ii)  if  J||  >  6. 

v)  Stop. 

Extensions  of  the  basic  problem  can  be  made  by  considering  the  following  decentral¬ 
ized  dynamic  stabilizing  compensator: 

it  —  -f  Gij/t 


Ui  =  H,Zi  +  M,yi 


i  =  1,. 


I/. 


The  case  of  dynamic  output  feedback  is  reduced  to  the  previous  problem  by  augment¬ 
ing  the  state  of  the  compensators  with  that  of  the  system,  i.e., 
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where  v,  =  ij  and  Ci  is  provided  by  the  designer;  modifying  the  cost  to 
J  f  +  2  X)  x'^NiUi  +  ^  ujR]ui  +  xjR]zi)  dt. 

^  «=i  «=i  «=i 

Further  extensions  of  the  above  problem  can  be  made  by  incorporating  frequency 
weighting  into  the  quadratic  criterion.  Centralized  frequency  weighting  has  been 
introduced  in  [5].  A  decentralized  frequency  weighted  cost  was  mentioned  in  [56]  and 
has  recently  been  further  developed  [59].  We  now  outline  some  of  the  steps  necessary 
in  reducing  the  decentralized  frequency  weighting  to  the  decentralized  LQR  problem. 

2.3.4  Decentralized  Frequency  Shaping  (DFS) 

Frequency  shaping  permits  frequency  dependent  weighting  matrices  in  the  quadratic 
cost  functional  [5].  This  allows  the  designer  to  use  some  of  the  classical  frequency 
domain  constraints  in  the  state-space  LQ  setting.  The  cost  functional  in  the  frequency 
domain  becomes: 

J  =  J  /_”  {x-Q^X  +  £  j  du,  (2.63) 

where  Qo  ^  0  and  Qi  and  Ri  for  e  =  l,...,i/  are  positive  semidefinite  and  positive 
definite  Hermitian  matrices  respectively. 

Frcqueticy  shaping  is  used  in  the  following  way.  As  presented  in  [60],  frequency  de- 
pendent  weighting  on  the  control  allows  for  better  robustness  properties  than  similar 
weighting  on  the  state,  so  the  control  input  is  chosen  for  the  frequency  dependent 
cost  weighting  and  the  output  weighting  remains  constant.  By  choosing  Ri{ju>)  to  be 
larger  for  large  uj  than  for  small  u;,  control  actions  are  penalized  at  high  frequencies 
more  than  low  frequencies.  Effectively,  this  choice  of  Ri{ji*})  reduces  the  loop  gain  of 
the  channel  at  high  frequencies  which,  in  turn,  results  in  a  system  that  can  better 
accommodate  high  frequency  uncertainties  related  to  channel  i. 

Using  the  notation  of  [59],  Qi  and  Ri  are  factored  as: 

Q,  =  Q’Qi 

Ri  =  k" Ri  i  =  I, ...  ,1/. 

For  this  decomposition,  the  cost  may  be  written  as 

^  [Y:Yi  +  (/•(/, ]^  dw  (2.66) 


(2.64) 

(2.65) 
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where 


Y,  =  Q,yu 

(2.67) 

Ui  =  RiUi  1  =  1,...,  I/. 

Assuming  the  Qi  are  proper  rational  matrices  and  the  Ri  are  rational  matrices  with 
entries  being  either  proper  or  non-causal,  they  can  be  realized  in  the  following  manner 
respectively: 

Vi  =  F'iVi  +  G^iVi 

(2.68) 

v;  =  Hfrfi  +  N^yi  *  =  l,...,i/ 

and 

i.  =  + 

b.  =  + 

p=l 

where  (•)^'’^  denotes  the  pth  derivative  of  (•)  with  respect  to  time.  Furthermore,  it 
is  assumed  that  the  entries  of  Ri{u})  have  k  more  zeros  than  poles.  For  additional 
results  in  decentralized  frequency  shaping,  see  [59]. 

2.3.5  Overlapping  Decompositions 

Due  to  the  complexity  and  large  scale  of  systems  describing  flexible  structures,  de¬ 
composition  methods  play  a  significant  role  for  problem  solving.  One  method  for 
controller  design  of  interconnected  s  ibsystems  is  a  decentralized  approach  which  as¬ 
sumes  that  each  subsystem  is  controlled  independently  using  local  information  and 
performance  considerations.  This  method  deals  with  disjoint  subsystems,  where  there 
is  no  information  overlap  in  either  the  dynamics,  control,  or  output.  However,  for  a 
large  class  of  systems  (e.g.,  flexible  structures,  traffic  regulation,  and  power  systems), 
the  sharing  of  information  among  controllers  is  essential.  For  systems  of  this  type,  it 
is  natural  to  consider  an  overlapping  decomposition  approach,  where  subsystems  can 
have  overlapping  of  states,  inputs,  and  outputs. 

The  underlying  concept  for  overlapping  decomposition  is  the  inclusion  principle.  In 
the  next  section,  we  introduce  the  inclusion  principle  and  several  concepts  relating 
to  inclusion  such  as  expansion/contraction  and  restriction/unrestriction.  Using  these 


(2.69) 
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roncepts,  controllers  are  designed  for  the  individual  (disjoint)  subsystems,  then  con¬ 
tracted  to  the  original  state-space.  Conditions  for  contraction  of  controllers  are  also 
discussed. 

The  Inclusion  Principle 

The  mathematical  framework  for  decentralized  control  using  overlapping  decompo¬ 
sition  is  the  Inclusion  Principle.  Here,  overlapping  subsystems  are  expanded  into  a 
larger  state  space  and  control  space,  such  that  the  subsystems  appear  to  be  decoupled 
in  state,  input,  and  output.  Control  laws  are  then  designed  for  these  smaller,  disjoint 
subsystems.  Using  certain  contraction  conditions,  the  subsystems  and  controllers  are 
contracted  to  the  original  space  for  implementation,  and  if  all  of  the  appropriate  con¬ 
ditions  are  met,  then  the  closed  loop  properties  that  were  derived  for  the  expanded 
decoupled  subsystems  are  intact  when  the  system  is  contracted.  For  this  structure, 
we  consider  the  state,  input,  and  output  inclusion  problem  for  controller  design. 

Consider  a  system  with  two  decentralized  control  agents  described  by 

X  =  Ax -i- B\Ui  +  B2U2,  i(0)  =  zo  6  IR" 

y,  =  Cix  (2.70) 

yi  =  C2X 

where  j  €  HI”  is  the  state,  Uj  €  IR”‘'  is  the  input,  and  y,  €  IR'’’  is  the  output  of  the 
zth  control  agent  {i  =  1,2).  It  is  assumed  that  the  system  given  above  does  not  have 
any  decentralized  fixed  modes  (DFM’s)  [61]. 

Consider  the  following  partition  of  the  state  x: 

^OlCj 

Xoc  (2.71) 

where  Xo,c,  €  IR”®*'’  is  observable  and  controllable  only  by  agent-i,  and  Xqc  €  IR”°‘ 
is  observable  and  controllable  by  both  agents.  Then  the  matrices  of  (2.70)  have  the 
following  forms  [62]: 

■  >1,,  0  O' 

A  =  A21  A22  ^23 

0  0  A33 
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A 

B  =  [B^  :  B2]  =  B21 

0 


(2.72) 


C=  ■ 


Cii  C12  0 

0  C22  C23 


Consider  the  transformation: 


X  =  Vx  = 


€  IR”,  n  =  n  +  rio 


(2.73) 


The  expansion  [Ij  of  the  original  system  (2.70)  with  respect  to  the  transforma¬ 
tion  (2.73)  is  given  by  (see  [9,  14]  for  developments): 


i  =  -I-  Bu, 

y  =  Cx 


where 


(2.74) 


x(0)  =  Vxo  €  /?" 


(2.75) 
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(2.76) 
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c  = 


Cn  Ci2 

0 

0 

A 

Cl 

C12 

0  0 

C22 

C23 

C21 

C2  . 

It  can  be  shown  [7]  that  the  response  x{t)  of  (2.70)  is  related  to  the  response  x{t) 
of  (2.74)  by 

a:(0  =  V^x{t)  V<  >  0  (2.77) 

where  is  the  generalized  inverse  of  V  satisfying  V^V  =  I.  Furthermore,  (2.70) 
and  (2.74)  have  equivalent  input/output  descriptions,  i.e., 

G{s)  =  CisI  -  A)-^B  =  G{s)  =  CisI  -  A)-'B  (2.78) 

The  systems  are  related  by  the  following  transformations  from  [9]: 

A  =  VAU  +  M 

B  =  VBQ  +  N  (2.79) 

C  =  TCU  +  L 

where  M,N  and  L  are  complimentary  matrices  of  appropriate  dimensions.  For  (2.74) 

to  be  an  expansion  of  (2.70),  a  proper  choice  of  M,  N  and  L  is  required  as  given  in  [9]. 

In  this  case,  the  transformation  matrices  are  given  as  follows: 
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Now  consider  the  design  model  described  by: 
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X  ■=  Ax  +  Bu 
y  =  Cx 


(2.80) 


where 


A 

B 

C 


Ai  0 
0  A2  > 

Bi  0 
0  B2 

Cl  0 
0  C2 


(2.81) 


The  transfer  function  matrix  (TFM)  for  the  design  model  (2.80)  is  given  by: 

G{s)  =  C{sl  -A)-^B  =  block  diag  [G,  (s),  G2(s)]  (2.82) 

where  G,(s)  is  the  TFM  for  the  decoupled  subsystem: 


Xi  —  AiXi  +  BiUi 

Vi  —  C,3;,-,  i  =  1,2 


(2.83) 


Note  that  (2.83)  is  the  minimal  local  model  for  agent-i. 


Model  Reduction 

The  decision  of  which  modes  to  retain  from  each  subsystem  can  be  made  by  con¬ 
sidering  the  cross  grammians  as  discussed  in  Section  2.2.2.  A  decentralized  control 
strategy  for  a  LSS  is  arrived  at  by  the  following  two  tasks:  (1)  decide  on  an  appro¬ 
priate  system  decomposition,  and  (2)  truncate  the  model  to  a  manageable  number  of 
modes.  Both  of  these  goals  can  be  realized  by  using  a  balanced  realization  [46].  In  the 
decentralized  case,  we  can  perform  a  decentralized  balanced  realization  for  each  I/O 
channel,  and  perform  modal  truncation  on  each  subsystem  by  deleting  those  modes 
that  are  weakly  controllable  (observable)  through  that  particular  I/O  channel.  Thus, 
a  decentralized  balanced  realization  not  only  gives  the  control  designer  insight  for 
model  reduction,  it  also  helps  to  decide  on  the  appropriate  subsystem  decomposition. 

In  [63],  Williams  presents  an  efficient  algorithm  for  computation  of  grammians  for 
models  given  in  modal  coordinates.  In  his  work,  decentralized  control  is  assumed, 
where  control  is  affected  through  a  given  I/O  pair  (uj,j/,),  or  u,  =  G,(s)y,(s). 


29 


2.3.6  Isolated  Boundary  Loading  and  Controlled  Component  Synthesis 
(CCS) 

Isolated  Boundary  Loading 

In  Component  Mode  Synthesis  (CMS)  [64,  65,  66],  a  large  flexible  structure  is  assumed 
to  be  composed  of  components  or  substructures.  Each  component  or  substructure  is 
separately  modeled  using  finite  element  modeling  techniques.  The  order  of  each  of 
these  component  models  is  reduced  using  the  assumed-modes  method  or  the  Rayleigh- 
Ritz  method  with  constraint  equations.  The  reduced  order  components  are  then 
mathematically  coupled  together  with  interface  compatibility  equations  to  form  a 
reduced  order  model  of  the  entire  structure. 

Component  modes  can  be  classified  cis  rigid-body  modes,  normal  modes,  constraint 
modes,  or  attachment  modes  [67,  68].  The  basic  idea  in  CMS  is  to  find  a  set  of  com¬ 
ponent  modes  for  each  component  that  make  the  greatest  contribution  to  the  modes 
of  the  composite  structure.  Several  methods  exist  to  select  these  modes.  Perhaps  the 
most  widely  known  methods  are  the  Craig-Bampton  method  [69],  the  MacNeal-Rubin 
method  [70,  71],  and  the  Benfield-Hruda  method  [72]. 

Isolated  boundary  loading  developed  by  Young  [13]  is  a  component  modeling  method 
which  is  based  on  the  boundary  loading  ideas  of  the  Benfield-Hruda  method.  Dynamic 
contributions  of  other  components  are  included  in  the  model  of  a  given  component. 

This  discussion  will  be  restricted  to  a  system  comprised  of  only  two  components. 
Stiperscripts  A  and  B  will  be  used  when  referring  to  quantities  associated  with  the 
two  components.  Each  of  these  components  is  £issumed  to  have  dimension  n.  In 
gciUTal,  any  number  of  components  containing  any  number  of  nodes  may  be  used. 

(.'M.S  methods  separate  the  component  nodes  into  two  groups,  interface  nodes  and 
internal  nodes.  Here,  the  internal  nodes  are  further  separated  into  internal  boundary 
nodes  and  internal  nodes.  The  internal  boundary  nodes  are  the  nodes  of  the  boundary 
element  that  are  not  actually  connected  to  the  adjoining  structure;  i.e.,  the  nodes 
that  connect  the  boundary  elements  to  the  internal  elements  of  the  same  component. 
The  remaining  nodes  are  the  internal  nodes.  Thus,  the  finite  element  model  of  a 
component  can  be  partitioned  as 

■  Mii  Mik  ®  1  [  9i  1  r  0  ■  r  g,  1  r  /,  1 

Mki  Mkk  Mkj  qk  +  AYi  f^kk  t^kj  qk  =  fk  (2 -04) 

_  0  Mjk  Mjj  J  [  gj  J  L  0  t<ik  J  L  J  L  fi  . 

where  g,,  g^,  and  g^  are  the  generalized  displacements  at  the  interface,  internal  bound¬ 
ary,  and  internal  nodes,  respectively. 
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In  isolated  boundary  loading,  a  truncation  matrix  is  used  to  eliminate  the  internal 
and  internal  boundary  nodes.  The  truncation  matrix  is  given  by 


^  nt 


0 

0 


(2.85) 


Thus,  the  interface  loading  matrices  are  merely  the  stiffness  and  mass  matrices  cor¬ 
responding  to  the  boundary  nodes  of  the  adjoining  component; 


bT"  i.'B 


A'"  =  A"  'I' 
A/®  =  . 


(2.86) 

(2.87) 


The  resulting  modified  stiffness  and  mass  matrices  for  component  A  are  given  by 
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By  letting  the  modified  boundary  submatrices  be 


i<L  =  K-Ui'H 


Equations  (2.90)  and  (2.91)  can  finally  be  written  as 


i<c,  A'4 

0  4\ 


(2.88) 

(2.89) 

(2.90) 

(2.91) 

(2.92) 

(2.93) 

(2.94) 
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(2.95) 


Similarly,  the  modified  stiffness  and  mass  matrices  for  component  B  can  be  written 
as 


Kt  = 
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where 
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(2.96) 

(2.97) 


(2.98) 

(2.99) 


Equations  (2.99)  and  (2.98)  indicate  that  the  modified  boundary  submatrices  are 
the  same  for  each  of  the  components  involved  and  may  be  found  either  by  loading 
component  A  onto  component  B  or  by  loading  component  B  onto  component  A. 


It  should  be  noted  that  the  modified  mass  and  stiffness  matrices  can  be  determined 
from  the  model  of  an  expanded  component  consisting  of  the  original  component  ex¬ 
tended  one  element  into  the  adjoining  component  [13].  The  model  of  the  expanded 
component  will  contain  the  internal  boundary  nodes  of  the  adjoining  component  in 
addition  to  all  of  the  nodes  of  the  original  component.  The  modified  mass  and  stiff¬ 
ness  matrices  are  then  found  by  deleting  the  rows  and  columns  associated  with  the 
internal  boundary  nodes  of  the  adjoining  component. 


The  coupling  of  substructures  using  isolated  boundary  loading  employs  superposition 
to  couple  the  components  at  the  boundary.  Essentially,  the  equation  of  motion  for 
the  boundary  of  the  coupled  structure  is  found  from  superposition  of  the  forces  and 
motions  at  the  boundaries  of  the  individual  components  A  and  B  such  that 


(2.100) 

(2.101) 

The  finite  element  model  for  the  coupled  structure  is  given  by 


=  qf  +  qf 


-f-  /C^  V®  =  Z'*® 


(2.102) 
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where 
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(2.103) 


(2.104) 


(2.105) 

(2.106) 


Using  isolated  boundary  loading  for  component  modeling  and  this  superposition- 
based  approach  for  substructure  coupling,  the  coupled  model  is  the  exact  model  for 
the  composite  system.  In  conventional  large-scale  system  design  methods,  the  interac¬ 
tions  between  components  are  generally  discarded  in  decoupled  substructure  models. 
The  component  models  formed  using  isolated  boundary  loading,  however,  include  the 
interactions  with  other  components.  The  component  models  are  identical  to  the  de¬ 
coupled  subsystem  models  obtained  using  overlapping  decomposition  (7,  8,  9,  10]  on 
the  composite  system  model.  Thus,  the  powerful  tools  used  in  overlapping  decompo¬ 
sition  analysis  of  large-scale  systems  may  be  applied  to  the  isolated  boundary  loading 
models. 


Controlled  Component  Synthesis 

Controlled  component  synthesis  (CCS)  is  a  control  design  method  developed  by 
Young  [13]  in  which  controllers  are  designed  for  individual  components  b^lsed  on  com¬ 
ponent  models  as  discussed  previously.  Using  conventional  large-scale  system  design 
methods  requires  a  full  model  of  the  entire  system.  This  model  is  typically  decom¬ 
posed  into  interconnected  subsystem  models  which  are  then  decoupled  by  ignoring 
subsystem  interaction  terms.  A  key  advantage  of  controlled  component  synthesis  is 
that  a  full  model  of  the  entire  structure  is  not  required.  This  decentralized  control 
strategy  extends  the  component  mode  synthesis  modeling  concepts  into  control  con¬ 
cepts  by  allowing  component  controllers  based  on  component  models  to  be  developed 
independently. 
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Overlapping  decomposition  can  be  used  to  derive  component  models  from  the  com¬ 
plete  model  of  a  structure.  Considering  the  two-component  system  example  from  the 
previous  section,  the  component  model  is  given  by 
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(2.107) 


The  second  term  on  the  right  side  of  this  equation  represents  the  dynamic  contribu¬ 
tions  of  the  internal  boundary  motion  of  component  B  to  the  motion  of  component  A. 
Neglecting  this  term  results  in  the  component  model  for  A  which  would  be  derived 
using  isolated  boundary  loading. 


The  interlocking  control  concept  is  one  of  the  key  ideas  behind  CCS.  With  inter¬ 
locking  control,  colocated  sensors  are  placed  at  the  internal  boundary  coordinates 
of  each  component.  The  component  controllers  are  designed  using  the  component 
models  as  derived  in  the  previous  section  to  minimize  motion  at  the  internal  bound¬ 
ary  nodes.  Essentially,  one  or  more  controllers  are  inserted  near  the  substructure 
interfaces  as  shown  in  Figure  2.1.  By  minimizing  the  internal  boundary  motion,  the 


Figure  2.1:  Inserting  a  Controller 

transmission  of  disturbances  between  components  is  suppressed,  the  second  term  in 
Efiuation  (2.107)  becomes  small,  and  the  accuracy  of  the  assumed  component  model 
infrea.ses. 
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Although  any  technique  may  be  used  for  controller  design,  the  linear  quadratic  op¬ 
timal  regulator  (LQR)  is  perhaps  the  most  convenient  [13].  The  weighted  internal 
boundary  states  are  to  be  included  in  the  cost  function  with  the  weighted  control 
inputs.  Because  the  sensors  are  located  at  the  internal  boundary  coordinates,  these 
states  are  merely  the  outputs  of  the  system.  The  control  inputs  are  u  =  /*  because 
the  actuators  are  also  located  at  the  internal  boundary  coordinates.  Using  the  LQR 
approach,  the  component  model  for  component  A  is  formed  as 
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(2.108) 


where  is  the  control  input.  Using  velocity  measurements,  the  output  is  defined  to 
be 


y^  =  9^  (2.109) 

The  control  law  for  component  A  is  found  to  minimize  the  cost  function 

J  dt.  (2.110) 


In  general,  a  state  feedback  control  law  could  be  found  to  minimize  this  performance 
criterion.  However,  optimal  output  feedback  can  also  be  employed  to  achieve  good 
results  and  is  easier  to  implement  than  state  feedback.  This  procedure  is  repeated 
for  each  component.  CCS  controller  designs  are  inherently  decentralized  since  only 
the  component  outputs  or  states  are  required  for  the  component  control  law. 


2.3.7  Maximum  Entropy/Optimal  Projection  (MEOP) 

One  of  the  fundamental  concerns  in  the  application  of  modern  control  designs  for 
active  vibration  suppression  is  the  loss  of  stability  and  the  performance  degradation  in 
face  of  modeling  errors.  This  is  especially  true  in  the  area  of  flexible  structures  which 
are  inherently  infinite-dimensional  eis  discussed  in  Section  2.1.1.  The  truncation  of  the 
infinite  dimensional  system  to  a  finite  dimensional  system  (as  in  the  finite  element 
representation  of  Section  2.1.2)  induces  truncation  errors  which  tend  to  be  more 
significant  at  high  frequencies.  Therefore,  the  theoretical  (calculated)  mode  shapes 
and  frequencies  are  less  accurate  for  higher-order  modes.  This  can  be  interpreted  as 
the  higher  modal  frequencies  and  mode  shapes  are  more  sensitive  to  small  details  of 
geometry,  construction,  and  material  properties.  The  high  frequency  modes  are  thus 
subjected  to  larger  errors  in  a  finite  element  model. 
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To  tackle  this  kind  of  control  problem,  Hyland  [73,  74,  75]  proposed  the  stochastic 
Maximum  Entropy  (ME)  approach  to  handle  model  uncertainties  in  order  to  improve 
the  system  robustness  with  respect  to  parameter  variations  and  high  order  unmodelled 
dynamics.  The  ME  approach  to  model  uncertainty  is  essentially  induced  by  the 
minimum  data  set  which  incorporates  as  little  statistical  parameter  information  as 
possible  consistent  with  the  available  a  priori  statistical  data.  The  basic  premise 
of  the  Maximum  Entropy  error  modeling  is  that  the  magnitude  of  the  error  is  a 
white  noise  process.  The  justification  for  the  use  of  the  Maximum  Entropy  stochastic 
approach  for  structural  control  design  is  reinforced  by  various  qualitative  results  for 
the  mean  square  optimization  problem  examined  in  [75]. 

From  the  view  point  of  control  engineers,  the  ME  approach  was  developed  as  a  means 
to  design  compensators  with  the  performance/robustness  trade-off  being  explicitly  de¬ 
termined  by  the  quadratic  cost  functional.  Thus,  the  ME  algorithm  is  essentially  a 
“robust”  extension  to  LQG  theory.  By  explicitly  including  the  parameter  uncertain¬ 
ties  in  the  model,  the  performance/robustness  trade-off  can  be  directly  manipulated 
to  suit  the  design  requirements. 

In  addition  to  the  potential  problem  of  modeling  errors  in  flexible  structures,  the  large 
dimensionality  of  the  dynamic  model  inherent  in  the  structural  system  poses  another 
problem  in  control  design  and  implementation.  It  is  known  that  the  classical  LQG 
controllers  have  the  same  order  as  the  controlled  plant.  If  the  plant  is  of  high  order, 
then  the  LQG  controller  dimension  may  be  too  large  for  real  time  implementation. 
Reduced  order  controller  design  has  therefore  received  much  attention  in  the  past  few 
years  for  flexible  structure  applications.  In  the  literature,  various  design  algorithms 
appear  which  involve  either  ad  hoc  reduction  and  truncation  schemes,  or  the  two-step 
approach  to  reduced  order  controller  design  in  which  a  compensator  order  reduction  is 
preceded  by  the  LQG  design  or  a  LQG  design  is  preceded  by  the  plant  order  reduction 
[46,  76,  77,  78,  79,  80].  However,  these  approaches  are  not  optimal,  and  worse  yet, 
they  may  cause  closed  loop  instability  [81].  The  Optimal  Projection  (OP)  approach 
proposed  by  Hyland  and  Bernstein  [82]  focuses  on  developing  direct  reduced  order 
design  methods  which  have  greater  flexibility  with  respect  to  the  dimension  of  the 
controller.  It  involves  the  minimization  of  a  quadratic  cost  functional  while  directly 
constraining  the  order  of  the  feedbcick  controller  to  guarantee  the  optimality  of  the 
reduced  order  design. 

In  order  to  design  a  controller  which  takes  care  of  both  the  robustness  issue  and 
the  controller  dimensionality  issue,  the  combined  Maximum  Entropy /Optimal  Pro¬ 
jection  (MEOP)  methodology  has  been  developed [83].  This  methodology  provides  a 
powerful,  yet  very  structured  algorithm  for  designing  a  robust,  fixed-order  dynamic 
compensator. 

Consider  an  order  linear  time  invariant  system 
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i  =  +  Bu  +  wi 


y  =  Cx  +  W2  , 


(2.111) 


under  the  usual  assumptions  of  controllability  and  observability,  where  W2  are 
zero  mean  white  noise  processes  with  noise  intensity  matrices  Vi  and  Vj,  respectively. 
It  is  required  to  design  an  order  robust,  zero  set-point  compensator 

Xc  =  AcX-\-Fy  u  =  —Kxc  (2.112) 

to  minimize  the  cost  functional 

J  =  Jim  {x^{t)Rix{t)  u^{t)R2u{t))dt^  (2.113) 

for  j:  G  IR",  €  IR’*%  u  €  IR”‘,  y  €  IR',  /I  €  1R"’‘”,  B  G  IR"’'”',  C  G  IR'^'", 
Ac  e  F  G  IR"'’^',  K  G  R,  G  1R"’‘"  (positive  semi-definite),  and 

R2  €  IR’"*"*  (positive  definite). 

The  first  order  necessary  conditions  for  the  quadratically  optimal,  steady  state,  ro¬ 
bust,  reduced-order  dynamic  compensation  are  the  existence  of  non-negative  definite 
matrices  P  G  IR"*",  Q  G  IR"*",  P  G  IR"*",  Q  G  IR"*"  satisfying  the  following 
coupled  Lyapunov  and  Riccati  equations: 

0  =  PAs  +  AlP  +  j2AjPAi  -  PJR:^}P,  +  Rx 

i=l 

+  DA  -  QXVC.fhAi  -  Q.Vf,'Q) 

i=l 

+t'[PBR2^B'^Pt^  (2.114) 

0  =  A,Q-\-  QAl  -  QsV2-'QI  +  V, 

1=1 

'ZiA,  -  B,R^^P,)Q{Ai  -  B,R-^lP,Y 

(2.115) 

0  =  PAqs  +  AlsP  +  PJR2}P,  -  tIPBR-'B'^Pt^  (2.116) 

0  =  ApsQ  +  QAls  +  Q,V2~'QI  -  rigC^Vj-'CQrJ  (2.117) 

where  the  proje''tion  operator  t  G  IR"*"  is  given  by 


T 


k=l  k 
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and  represents  the  A;"*  eigenprojection  of  QP,  with  =  In  —  t.  In  (2.114)- 

(2.117),  Ai  6  Bi  €  Ci  €  IR^^”  are  the  uncertainty  matrices  of  the  plant 

[S3],  from  which 

As  =  A  +  \Et=i 

Cs^C-¥\Y,UCiAi 

V2,  =  V^2  +  Er=i  Ci{Q  +  Q)CJ 
Ps  =  Bjp  +  Er=i  BJ{P  +  P)Ai 
Qs  =  QCJ  +  Er=i  A,{Q  +  Q)CT 


where  /z  is  the  number  of  sets  of  uncorrelated  uncertainties. 

In  terms  of  the  solution  to  the  above  MEOP  equations,  the  compensator  dynamics 
are  specified  according  to 


(2.118) 

F  =  VQaVia^ 

(2.119) 

K  = 

(2.120) 

where  the  operators  G  €  IR”‘‘^'‘,r  G  must  satisfy  = 

and  G^r  =  t. 

2.3.8  Hoc  Design  Approach 

A  major  emphasis  in  control  design  is  the  issue  of  robustness.  As  we  have  men¬ 
tioned  previously,  uncertainties  tend  to  exist  in  system  models.  In  the  case  of  a  large 
flexible  system,  there  are  uncertainties  in  parameters  such  ais  dampings  and  frequen¬ 
cies  of  flexible  modes.  Also,  models  of  flexible  structures  typically  neglect  dynamics 
of  higher  frequency  modes.  Our  control  objectives  must  be  achieved  despite  these 
modeling  errors.  These  two  types  of  uncertainties  can  be  described  as  structured 
and  unstructured  uncertainty.  The  control  system  must  be  robust  enough  to  handle 
both  types  of  uncertainty  and  still  be  able  to  keep  the  designed  closed-loop  perfor¬ 
mance.  The  control  method  described  here  is  an  //cx>-baised  technique  of  “weighted 
gap  optimization,”  which  is  closely  related  to  Glover-McFarlane  loop-shaping.' 


5.  =  B+|Er=i  AiBi 

R2s  =  R2-\-  y:u  bJ{p  -b  p)B. 

Aqs  =  a,  — 

Aps  =  A,  —  BaR^}  Pa 


Weighted  Gap  Optimization 

The  (weighted)  gap  metric  is  a  measure  of  uncertainty  b£ised  on  the  notion  of  the 
graph  of  a  dynamical  system.  The  primary  motivation  for  this  distance  measure  is 

‘The  Hro  work  presented  here  was  performed  by  Scott  Buddie  for  his  M.S.  thesis  [84]  under  the 
guidance  of  Professor  Malcolm  Smith  and  was  not  directly  supported  by  this  contract.  The  work 
i.s  inrhidefl  for  completeness. 


38 


Figure  2.2:  Standard  Feedback  Configuration 


the  need  to  handle  unstable  and  marginally  stable  systems  in  a  satisfactory  way  from 
the  point  of  view  of  uncertainty.  It  may  be  claimed  that  a  graph-based  distance 
measure  has  certain  advantages  in  other  situations  also,  such  as  when  lightly  damped 
modes  are  present,  as  is  the  case  with  flexible  structures. 

The  gap  between  systems  P  and  Pi  can  be  defined  as 


^(P,P,):=max{^  (P,P,),6  (Pi,P)} 

where  the  directed  gap  S  (P, Pi)  is  given  by 


l{P,Px) 


(2.121) 


(2.122) 


and  P  =  NM~^  and  Pi  =  N\Mi^  are  normalized  right  coprime  factorizations  of 
the  system  transfer  functions  [85].  For  weighting  functions  Wi,Wo  €  Hoo  such  that 
VF"’,  €  Hoo,  the  corresponding  weighted  gap  metric  6{P,  Pi,Wo,Wi)  is  defined 

as: 


6{P,Pi,Wo,Wi)  :=  6{WoPWi,WoPiWi).  (2.123) 

By  using  these  weighting  functions,  one  can  emphasize  or  de-emphasize  differences 
between  the  two  plants  at  desired  frequencies  (86). 

Consider  the  feedback  configuration  of  Figure  2.2  where  P  and  C  are  linear  systems. 
The  configuration,  denoted  by  [P,C]  is  defined  to  be  stable  if  the  operators  Xi  —*  Cj 
for  i,j  =  1,2  are  bounded.  This  is  equivalent  to 
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(2.124) 


(I  -c]  f  (I-  CP)-^  C(1  -  PC)-^  \ 

\-p  I  )  ■  V  P(/  -  CP)-^  (/  -  PC)-'  ) 

being  in  Poo- 

Let  [P,  C]  be  a  stable  feedback  configuration  and  define 

KP,  c,  ►*'i)  •■=  U  JJTp  )  (/  -  CP)-'  (w.  -  cw:') 

as  the  maximal  amount  of  uncertainty  in  the  weighted  gap  that  a  feedback  system 
can  tolerate  [87,  88]. 

The  optimal  robustness  radius  can  be  defined  as 

6opt(P)  :=  max6(P,C,W^o,W^.). 

C  stbl2 

or 

(ci^K  zCp  ) 

The  problem  of  weighted  gap  optimization  reduces  to  solving  this  well-posed  Hoo- 
minimization  problem  [85,  89]. 

The  problem  of  finding  a  controller  that  is  optimally  robust  with  respect  to  the 
gap  metric  is  identical  to  the  problem  concerning  robust  stabilization  under  coprime 
fraction  uncertainty.  This  implies  that  any  technique  used  to  find  a  controller  for 
coprime  fraction  uncertainty  can  be  used  to  solve  the  gap  metric  problem.  A  solution 
to  the  coprime  fraction  uncertainty  problem  is  given  by  Glover  and  McFarlane  [90, 
91,  92]. 

In  the  work  of  McFarlane  and  Glover  [92],  a  certain  “loop  shaping”  design  proce¬ 
dure  is  proposed.  This  procedure  begins  with  an  augmented  plant  transfer  function 
Paug  •=  WoP^i  Containing  input  and  output  “weighting”  functions  Wg  and  VFi.  The 
weighting  functions  are  selected  on  the  basis  of  a  desired  loop  shape  for  the  control 
system.  The  controller  is  designed  to  give  optimal  robustness  with  respect  to  normal¬ 
ized  coprime  fraction  perturbations  of  Paug^  We  point  out  that  this  choice  of  design 
synthesis  problem  is  precisely  the  same  as  the  .^oo-minimization  problem  given  by 
Equation  (2.126).  This  algorithm  was  primarily  used  in  calculating  the  controller  for 
the  truss  structure. 

Safonov  and  Chiang  claim  that  Moore’s  reduced  model  can  be  constructed  without 
computing  a  balancing  transformation.  This  can  be  done  by  constructing  projections 


-1 

(2.125) 

OO 
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in  terms  of  bases  for  the  left  and  right  eigenspaces  of  the  product  of  the  controllability 
and  observability  grammians  EcSo  [93].  There  are  two  methods  which  may  be  used  to 
perform  the  first  step  of  this  algorithm  which  computes  matrices  whose  columns  form 
bases  for  the  left  and  right  eigenspaces  of  EcEo  associated  with  the  “big”  eigenvalues. 
The  first  method  employs  the  Schur  form  of  the  matrix  Ec^o  and  a  Givens  rotation 
to  produce  ordered  values  along  the  diagonal  of  the  matrix.  The  second  method  is 
based  on  using  the  “square-roots”  of  Ec  and  Eo.  We  chose  to  use  the  first  method. 
The  decision  to  truncate  modes  is  made  by  choosing  the  desired  number  of  columns 
of  the  matrices  whose  columns  form  the  bases  for  the  left  and  right  eigenspaces.  It  is 
important  to  note  that  the  resulting  transformed  system  is  not  a  balanced  realization, 
but  it  is  an  equivalent  realization  of  the  balanced  realization  obtained  by  Moore’s 
algorithm.  Therefore,  if  a  system  needs  to  be  balanced,  this  approach  cannot  be 
used.  Our  concern  in  this  effort  was  to  avoid  numerical  conditioning  problems,  so  an 
equivalent  realization  was  deemed  sufficient. 
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3.  EXPERIMENTAL  HARDWARE  CONFIGURATION 


3.1  Truss  Structure,  Controllers  and  Actuators 

The  structure  is  a  12-meter  truss  [43,  44]  oriented  in  the  vertical  cantilever  position, 
as  shown  in  Figures  3.1  and  3.2.  It  is  composed  of  four  sections  with  4  truss  bays 
each.  These  four  sections  are  bolted  together  to  form  the  truss  with  a  total  of  16 
bays.  The  truss  is  constructed  from  square  aluminum  tubing,  with  cross  members  of 
Lexan  plastic  bolted  into  place. 

The  active  damping  is  achieved  through  eight  momentum  exchange  proof-mass  ac¬ 
tuators.  These  devices  are  essentially  linear  DC  motors  with  permanent  magnets 
mounted  on  low  friction  sliding  shafts  and  coils  fixed  to  the  housings  as  shown  in  Fig¬ 
ures  3.3  and  3.4.  The  proof  mass  is  centered  with  lightweight  springs.  Each  actuator 
is  supplied  by  its  own  current  driver,  which  also  incorporates  viscous  damping  via  a 
linear  velocity  transducer  (LVT)  mounted  between  the  proof  mziss  and  the  actuator 
housing.  The  actuators  are  capable  of  supplying  up  to  1  pound  force  from  2  Hz  to 
100  Hz.  The  locations  of  the  actuators  are  shown  in  Figure  3.5.  There  are  pairs  of 
actuators  mounted  perpendicular  to  each  other  at  the  center  of  both  the  one-half  and 
three-quarter  stations  on  the  truss  as  shown  in  Figure  3.6.  These  actuators  affect  only 
the  bending  modes,  not  the  torsional  modes.  There  are  four  actuators  at  the  tip,  two 
in  each  direction,  mounted  off  center  as  shown  in  Figuretipstat.  This  orientation  is 
used  to  affect  both  torsion  and  bending.  Actuators  can  be  monitored  via  three  video 
cameras  located  on  scaffolding  beside  the  truss.  The  scaffolding  also  provides  easy 
access  to  any  of  the  actuator  stations.  The  overall  system  is  shown  in  Figure  3.8. 
Sensors  for  the  control  system  are  accelerometers,  which  are  mounted  on  the  actuator 
housings.  The  acceleration  signals  are  integrated  using  bi-quad  integrators  with  a 
break  frequency  of  0.1  Hz  to  provide  velocity  signals.  A  point  light  source  is  located 
at  the  tip — 24  inches  from  the  center  of  the  truss.  A  digital  camera  can  record  the 
position  of  the  light  source,  and  thus  provide  information  regarding  the  displacement 
of  the  tip. 

During  testing,  the  truss  is  excited  via  two  alternative  sources.  The  first  alternative 
is  an  additional  proof-mass  actuator,  identical  to  the  control  actuators.  This  adds 
dynamics  to  the  structure  which  have  been  taken  into  account  in  the  model  of  the 
structure.  The  second  alternative  is  a  large  shaker,  mounted  to  the  wall  at  the  tip 
of  the  truss.  The  shaker  applies  a  force  to  the  truss  through  a  light  connecting  arm, 
and  is  capable  of  delivering  forces  up  to  10  lb. 
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Figure  3.1:  A  12-Meter  Vertical  Cantilever  Truss  Structur 


Figure  3.2:  12-Meter  Truss  Experiment  in  Vibration  Test  Facility 
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Figure  3.3:  Momentum  Exchange  Actuator 
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Figure  3.4:  12-Meter  Truss  Actuator 


nr  STATION  ACTUATOR  CONnaUlWnON 


Figure  3.5:  Actuator  Locations  on  Truss 
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Figure  3.7:  12-Meter  Truss  Tip  Station  Actuators 


Figure  3.8:  System  Overview 


3.2  Real-time  Control  System 


The  control  algorithms  are  implemented  digitally  using  a  custom  Optima  S  control 
computer  manufactured  by  Systolic  Systems,  Inc.  The  system  consists  of  the  control 
computer  and  a  Sun  3/50  graphics  workstation  that  serves  as  a  development  system. 
The  control  computer  has  a  Motorola  68030  host  processor  with  a  25  MHz  clock 
and  4  MB  of  RAM.  For  high  speed  vector  calculations,  the  Optima  S  has  a  systolic 
array  processor  with  a  peak  rating  of  32  million  floating  point  operations  per  second. 
External  control  is  accomplished  through  a  12-channel,  16-bit,  A/D  and  D/A  card. 
The  control  programs  are  written  in  C  programming  language,  compiled  on  the  Sun 
workstation,  and  downloaded  to  the  control  computer.  Unfortunately,  the  interface 
between  the  Optima  3  and  Sun  is  extremely  slow. 

Nevertheless,  this  system  is  capable  of  implementing  rather  complex  control  algo¬ 
rithms  at  sufficient  frequencies  for  truss  control.  For  example,  a  14-state  dynamic 
controller  can  easily  be  run  at  sampling  periods  less  than  5  ms. 

Data  acquisition  and  analysis  can  be  performed  either  by  uploading  data  to  the  Sun 
host  or  through  an  OnoSokki  dual  channel  spectrum  analyzer.  Since  the  interface 
between  the  control  computer  and  the  workstation  is  very  slow,  the  spectrum  analyzer 
is  frequently  used  to  minimize  the  time  required  to  collect  data  for  each  experiment. 
Usually,  rms  values  of  the  sensor  and  actuator  signals  were  calculated  on  the  Optima 
3  to  take  advantage  of  the  array  processing  capabilities.  The  rms  values  were  used  as 
a  quick  performance  summary  of  a  given  experiment  for  determining  whether  or  not 
to  save  (upload)  the  entire  data  record. 
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4.  SOFTWARE  FOR  CONTROL,  SIMULATION,  AND  DATA 

ACQUISITION 


4.1  Goals  of  the  Software 

One  of  the  aims  of  the  AF  in-house  project  was  to  develop  a  real  time  computer 
control  system  which  could  be  programmed  in  a  higher  level  language.  This  would 
provide  several  benefits: 

•  Real-time  control  tasks  such  as  reading  A/Ds,  controlling  sampling  period,  etc. 
are  all  performed  in  easy  to  access,  prewritten  subroutines. 

•  Controller  implementation  is  made  easier — control  laws  can  be  performed  in 
matrix  form. 

•  Ability  to  store  data  during  experiment,  then  upload  for  post  processing. 

•  User  friendly  software  possible. 

As  a  result,  the  control  and  simulation  software  written  can  be  used  by  an  individual 
with  minimal  experience  with  the  computer  system.  By  the  same  token,  additional 
controllers  could  be  added  to  the  software  with  minimal  understanding  of  C  program¬ 
ming  language  and  specific  knowledge  of  the  Sun  and  Optima  3  computers. 


4.2  Control/Simulation  Program 

Both  control  and  simulation  algorithms,  including  several  different  types  of  con¬ 
trollers,  are  implemented  entirely  in  one  piece  of  software.  This  was  done  for  a 
number  of  reasons.  First,  we  wanted  to  be  able  to  test  the  controller  software  as 
thoroughly  as  possible  without  actually  controlling  the  truss  or  running  the  software 
on  the  Optima  3.  This  enables  the  user  to  write,  compile,  and  test  the  software  on 
any  Sun  workstation.  Second,  it  was  felt  that  it  would  save  effort  in  the  end,  since 
changes  to  the  controllers  or  program  features  would  only  have  to  be  made  to  one 
program.  Third,  it  also  gives  the  benefit  of  being  able  to  run  simulations  on  the 
Optima  3,  which  greatly  decreases  the  run  time. 

4.2.1  Software  Overview 

The  control/simulation  software  is  interactive — prompting  the  user  for  various  pa¬ 
rameters  about  the  run.  A  sample  run  of  the  program  is  given  in  Appendix  A.  The 
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prompts  should  be  self  explanatory.  There  has  been  some  effort  to  make  the  program 
accept  only  correct  responses,  although  it  is  certainly  possible  to  enter  erroneous  data. 
Throughout  the  run,  the  program  reports  various  events  to  the  user  each  preceded  by 
an  ellipsis  (...).  These  reports  are  merely  intended  as  reminders  of  what  is  happening 
and  as  an  aid  for  debugging. 

The  software  has  been  divided  into  five  files  to  make  editing  a  bit  more  manageable. 


exp.c  This  is  the  main  program,  which  contains  all  the  prompts  to  the  user  and 
performs  all  the  calculations  for  control. 

subroutines.c  This  file  contains  several  important  subroutines.  The  getjuatrix 
routine  loads  ASCII  files  containing  floating  point  numbers  into  array  variables. 
The  fflvffl  routine  multiplies  a  matrix  and  a  vector.  This  routine  uses  the  array 
processor  when  compiled  for  the  Optima  3  and  is  one  of  the  most  important 
routines  for  implementing  controllers. 

band-vhite.iir  .c  This  routine  generates  a  bandwidth  limited  white  noise  signal. 
The  routine  uses  the  r2uid  command  to  generate  random  short  integers  which 
are  summed  to  produce  a  Gaussian  distribution.  The  resulting  signal  is  then 
filtered  with  a  Finite  Impulse  Response  (FIR)  filter  to  limit  the  bandwidth  to  50 
Hz.  If  another  cutoff  is  desired,  another  FIR  filter  could  easily  be  designed  using 
MATLAB.  The  filter  changes  the  RMS  value  of  the  signal.  To  compensate  for 
this,  the  standard  deviation  of  the  signal  is  calculated  after  filtering  and  then 
corrected  to  the  desired  value. 

sinjdisturbance .  c  This  routine  generates  a  sinusoidal  disturbance  and  has  the  same 
form  as  band-vhite-fir.c 

sin^veep.c  This  routine  generates  a  sinusoidal  disturbance  with  a  frequency  which 
“sweeps”  over  some  user-specified  frequency  range.  This  routine  is  similar  to 
sinjdisturbance. c. 

There  is  one  more  file,  dumoyjnain.c,  which  is  used  to  get  around  a  programming 
problem  with  the  A/D  and  D/A  routines  provided  by  Systolic.  The  A/D  and  D/A 
routines  apparently  need  to  be  aligned  on  16-bit  memory  words  instead  of  8-bit  mem¬ 
ory  words;  if  they  are  placed  after  other  routines  in  memory,  the  compiler  may  not 
align  them  correctly.  In  C,  the  main  routine  must  appear  first.  To  get  around  this, 
a  false  or  “dummy”  main  is  used  in  dummy jiain.c.  This  false  main  appears  first  in 
memory,  but  merely  calls  the  real  main,  contained  in  exp.c,  which  appears  after  the 
A/D  and  D/A  routines. 
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Although  the  software  may  appear  to  be  long,  it  is  relatively  straightforward.  Com¬ 
ments  appear  throughout  the  code.  The  major  tasks  performed  by  the  software  are, 

1.  Appropriate  matrices  are  loaded  from  data  files 

2.  Memory  is  allocated  for  signal  storage 

3.  Disturbance  signal  is  pre-calculated 

4.  Sensor  biases  are  checked  (if  controlling  truss) 

5.  Control/Simulation  loop  run 

6.  RMS  signal  values  calculated 

7.  .m  file  written,  containing  report  on  the  run 

8.  .dat  file  written,  containing  various  sampled  data  signals  from  truss 


4.2.2  Computer  Compatibility 

The  program  is  written  in  C,  using  several  standard  libraries  and  several  special  ones. 
For  most  common  operations,  the  standard  std.l,  stdio.h,  and  math.h  libraries  are 
used.  For  matrix  operations,  the  Mathpac  routines  provided  by  Systolic  systems  are 
used.  Systolic  has  also  provided  routines  to  read  from  the  A/Ds,  write  to  the  D/ As, 
control  the  sampling  period,  and  measure  intervals  of  time.  These  routines  are  in  an 
object  file,  cvrtl6.o. 

The  program  is  compiled  on  a  Sun,  but  it  can  be  compiled  to  run  on  three  different 
types  of  machines,  using  the  following  makefiles, 

•  Sun  3/50  or  3/60  —  makefile 

•  Sun  3/260  —  make^honon 

•  Optima  3  —  make-optima 


Different  versions  of  the  std.l,  stdio.h,  and  math.h  libraries  must  be  used  when 
compiling  for  the  Optima,  and  these  libraries  have  been  provided  by  Systolic.  The 
Mathpac  routines  are  also  different,  since  they  are  written  to  run  on  the  array  pro¬ 
cessor.  Only  a  small  number  of  Mathpac  routines  arc  available  to  run  on  the  array 
processor,  namely  vadd,  vsub,  vmul,  vdiv,  vsadd,  vsmul,  vmsa,  and  vsabs.  These 
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routines  are  in  the  library  qc.l.  It  would  be  possible  to  run  all  of  the  Mathpac 
routines  on  the  Optima  3’s  host  processor  at  a  reduced  speed  instead  of  the  array 
processor  if  the  Mathpac  source  code  was  available. 

Compiler  directives  are  used  to  make  the  program  compile  with  the  correct  libraries, 
calls  to  subroutines,  etc.  For  example,  if  the  program  is  not  compiled  to  run  on  the 
Optima  3,  the  A/D  and  D/A  routines  are  not  included,  since  the  Sun  will  not  have 
this  capability.  When  changing  from  one  machine  to  another,  all  of  the  program 
object  files  must  be  erased.  This  will  insure  that  all  the  files  are  recompiled  to  run 
on  the  correct  system.  If  this  is  not  done,  the  linker  will  not  be  able  to  find  all  the 
subroutines  called. 

4.2.3  Memory  Usage  and  Allocation 

The  Optima  3  has  two  types  of  memory.  There  is  4  MB  of  main  memory,  which 
can  be  accessed  only  through  the  host  processor.  This  memory  is  used  to  store 
the  program  and  most  of  the  variables.  When  running  a  C  program  on  the  host, 
variables  are  automatically  placed  in  this  memory  as  well  as  locations  reserved  with 
malloc  and  calloc.  Unfortunately,  there  does  not  appear  to  be  any  way  to  check 
how  much  memory  is  used.  If  too  much  is  requested,  the  program  craishes.  The 
executable  program  requires  about  160  KB  and  memory  is  allocated  for  variables  as 
it  is  requested  by  the  user.  This  limits  the  amount  of  memory  the  program  actually 
uses.  Most  of  the  memory  is  used  to  store  the  real-time  signals  measured  during  the 
experiment.  The  amount  of  memory  available  for  storing  real-time  signals  from  the 
experiment  allows  at  least  30  seconds  of  data  at  10ms  sampling.  The  usefulness  of 
this  storage  is  limited,  however,  since  the  data  link  between  the  Optima  3  and  Sun 
host  is  slow. 

The  second  type  of  memory  is  on  the  array  processor  card.  It  is  the  only  memory 
that  the  array  processor  can  read,  although  the  host  can  read  this  memory  also.  This 
memory  must  be  allocated  with  the  APalloc  command,  which  returns  a  pointer  to  a 
floating  point  number.  All  variables  used  in  calculations  on  the  array  processor  must 
be  stored  in  variables  allocated  in  this  manner.  Only  32-bit,  floating  point  numbers 
may  be  used  with  the  array  processor. 

4.2.4  Sampling  Rate 

The  sampling  rate  of  the  controller  is  precisely  controlled  using  the  A/D  and  D/A 
routines  provided  by  Systolic.  The  sampling  rate  is  set  when  the  initialization  routine 
CVTJNIT  is  run.  The  routine  CVT_STRT  starts  the  sampling  rate  timer.  The  first  A/D 
call  and  the  first  D/A  call  each  start  separate  timers.  When  the  software  reaches 
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Figure  4.1 :  Speed  of  a  Matrix/Vector  Multiplication  using  Optima  3’s  Array  Processor 

the  second  A/D  and  D/A  calls,  the  program  waits  until  the  proper  time  has  elapsed 
since  the  previous  calls.  The  program  will  stop  and  print  an  appropriate  message  if 
this  sampling  rate  is  exceeded.  There  is  a  two  sampling  period  delay  in  the  whole 
closed-loop  cycle. 

When  compiled  to  run  on  the  Optima  3,  the  program  makes  use  of  the  array  processor 
for  matrix/vector  multiplications  and  vector  additions.  For  the  multiplication,  the 
Mathpac  routine  vmsa  is  used.  It  performs  the  dot  product  of  two  vectors.  Due  to 
the  pipeline  architecture  of  the  array  processor,  this  operation  can  be  performed  at 
virtually  the  same  speed,  regardless  of  the  size  of  the  vectors.  This  routine  is  called 
once  for  every  row  in  the  matrix,  thus  for  square  matrices,  the  speed  of  execution 
will  be  proportional  to  the  order,  rather  that  the  order  squared.  Figure  4.1  shows  the 
level  of  performance  achieved. 

It  is  easy  to  see  that  the  matrix  equations  should  be  written  with  the  smallest  possible 
number  of  rows  in  the  equations.  Diagonal  matrices  can  be  multiplied  with  one  vmsa. 

The  program  as  written  can  sample  the  A/D,  write  to  the  D/A,  store  all  the  signals 
and  take  care  of  all  the  other  overhead  in  about  1 .5  ms.  A  14-state  dynamic  controller 
can  be  implemented  with  a  sampling  period  of  5  ms.  Perhaps  some  work  could  be 
done  to  reduce  this.  • 
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4.2.5  Output  of  the  Software 


The  time  response  data  generated  in  the  program  is  written  to  a  data  file  in  buffered 
binary  output,  which  can  be  read  directly  by  MATLAB.  A  MATLAB  header  is  written 
to  the  file  before  any  data  to  configure  the  file  as  a  MATLAB  .mat  file  [94].  The 
rate  at  which  data  points  are  written  is  selected  by  the  user,  thus  the  file  size  can  be 
controlled.  Data  are  written  to  the  file  one  data  block  at  a  time  using  the  C  command 
fwrite.  Using  the  fwrite  command  decreases  the  time  required  to  upload  data  to 
approximately  one  third  to  one  fifth  the  time  required  using  fprintf .  The  size  of  the 
data  blocks  vary  depending  on  the  machine.  The  data  written  to  the  file  may  consist 
of  any  of  the  following:  the  time  of  the  sampling  instant,  eight  sensor  outputs,  eight 
controller  outputs,  nine  mode  displacements,  force  to  shaker,  disturbance  signal  and 
the  LOS  pointing  error.  Usually,  some  subset  of  these  signals  are  written  to  the  file 
because  of  the  time  required  to  upload  the  data.  This  file  is  placed  in  a  subdirectory 
called  ./data. 

In  addition  to  the  data  file,  a  MATLAB  command  (.m  file)  file  is  written,  containing 
information  about  the  run.  This  file,  in  addition  to  the  data  file,  is  intended  to  be 
used  with  the  post  processor,  which  is  described  in  Section  4.3.  However,  the  .m  also 
gives  a  good  report  of  the  rms  values  of  the  sensor  signals,  which  can  be  printed  as  a 
record  of  the  run  or  used  for  a  quick  analysis  of  the  results. 

4.2.6  Addition  of  New  Controllers 

If  additional  controllers  are  designed  for  the  truss,  whether  by  OSU  or  another  party, 
there  is  no  need  to  write  additional  software.  New  controllers  can  be  added  to  the 
existing  software  relatively  easily.  Minimal  knowledge  of  C  programming  is  needed.  A 
good  reference  is  [95],  since  it  was  used  by  Systolic  Systems  Inc.  when  they  wrote  the 
C  language  routines  for  the  Optima  3.  It  is  also  very  helpful  to  read  the  documentation 
files  written  by  Systolic,  found  in  /usr/download. 

The  controller  should  be  written  in  matrix  form,  to  enable  use  of  the  mvm  and  vadd 
routines  which  run  on  the  array  processor.  The  constant  matrices  should  be  written 
in  ASCII  files,  one  line  per  row  (for  readability  only).  Any  format  readable  by  the  Xf 
option  of  the  C  command  f  scanf  is  acceptable,  but  it  is  advisable  to  use  exponential 
format  and  «is  many  decimal  places  as  possible.  Each  matrix  should  be  in  a  separate 
file.  These  files  can  be  generated  directly  from  CtrlC  (using  the  command  save  var 
>filen8une  -r  -c  -f  -(lx,  Ip5el5.7)  )  or  MATLAB  (using  the  command  save, 
then  tremslate).  If  the  controller  is  dynamic,  it  must  be  discretized  for  the  sampling 
period  at  which  it  is  to  be  run  before  it  is  saved  in  these  files. 

The  program  needs  modifications  in  the  file  exp .  c  only.  It  is  easy  to  spot  the  necessary 
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modifications,  by  studying  the  implementation  of  various  controllers  already  coded. 
Basically,  there  are  five  areas  to  modify: 

1.  Define  pointers  to  new  variables  used  in  the  control  law,  making  sure  they  are 
not  already  defined.  The  variables  used  may  be  the  same  £is  in  other  control 
laws,  since  only  one  controller  is  run  at  a  time,  however,  they  cannot  be  defined 
twice. 

2.  Add  controller  to  main  menu,  giving  it  a  controller  number,  which  will  be  its 
reference  throughout  the  program. 

3.  Add  a  case  to  the  switch  statement  that  loads  the  controller  files.  Also,  at  this 
point,  memory  for  the  constant  matrices,  controller  states,  and  any  temporary 
variables  is  allocated.  If  the  controller  is  dynamic,  there  may  be  a  switch  to 
load  controllers  discretized  at  different  sampling  rates.  The  sign  of  any  matrices 
should  be  corrected  so  there  are  no  subtractions  in  the  matrix  equations. 

4.  Similarly,  a  switch  must  be  added  to  calculate  the  feedback  signal  in  the  control 
loop.  This  relates  the  plant  output  vector  y  to  the  control  vector  u.  If  possible, 
use  the  mvm  and  vadd  routines  only.  Since  these  only  have  two  operands  and  they 
cannot  be  nested,  temporary  variables  may  need  to  be  used.  Also  remember  that 
all  of  these  variables  must  be  allocated  on  the  array  processor  using  APalloc. 

5.  Add  a  case  in  the  switch  that  prints  the  controller  title  to  the  .m  file. 

4.3  Post-Processor  Routines 

The  post-processor  is  a  menu  driven  plotting  routine  written  especially  for  the  sim¬ 
ulation  program  in  MATLAB  command  files.  It  is  intended  that  anyone  with  even 
minimal  knowledge  about  MATLAB  be  able  to  run  the  post-processor.  The  post¬ 
processor  can  bring  information  about  the  simulation  from  the  C  program  and  write 
it  on  the  plots,  eliminating  possible  mix  up  of  files. 

Unfortunately,  the  low  speed  of  the  interface  between  the  Optima  3  and  the  Sun 
has  severely  limited  the  usefulness  of  this  program.  With  a  feister  interface  it  would 
be  more  useful.  In  any  case,  the  following  sections  provide  a  short  summary  of  its 
capabilities. 


4.3.1  Capabilities  and  Structure 

The  post-proces.sor  is  menu  driven,  calling  many  individual  subroutines  for  specific 
task.s.  It  can  graph  actuator  signals,  sensor  signals,  fast  Fourier  transforms  of  sensor 
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signals,  pointing  error  and  mode  displacements.  The  user  can  view  these  graphs  with 
an  option  to  print,  or  print  them  out  directly.  Menu  options  exist  to  plot  groups  of 
graphs.  Additional  features  could  be  added  easily.  A  sample  run  of  the  post-processor 
is  shown  in  Figure  4.2. 


>>  Sim.l 

»  pp 

Post-Processor  for  12  m  Truss  Simulations 

1)  View  /  Plot  Actuator  Signals 

2)  View  /  Plot  Sensor  Signals 

3)  View  /  Plot  Actuator  Power 

4)  View  /  Plot  FFT  of  Sensor  Signals 

5)  View  /  Plot  Mode  Displacements 

6)  View  /  Plot  RMS  Pointing  Error 

7)  Plot  All  Actuator  Signals  and  Total  Power 

8)  Plot  All  Sensor  Signals 

9)  Plot  All  Actuator  Signals  emd  Sensor  Signals 

10)  Plot  All  Mode  Displacements 

11)  Plot  Total  Power  and  RMS  Pointing  Error 
99)  Exit  to  matlab  -  Load  new  data 

Option  ?  1 

Enter  Actuator  #  or  a  for  all  Actuators  1 
. . .  Plotting 

Do  you  weint  hardcopy  ?  (y/n)  n 

Figure  4.2:  Sample  Run  of  Post- Processor 

Each  graph  is  labeled  with  the  parameters  that  are  passed  from  the  simulation  pro¬ 
gram.  This  is  done  by  writing  a  file  that  will  be  executed  by  MATLAB  as  a  command 
(.m  file).  This  file  must  be  executed  prior  to  running  the  post-processor  routines,  as 
it  also  runs  an  initialization  routine. 

The  program  was  originally  written  to  plot  all  the  sensor  and  actuator  signals.  It 
is  clear,  however,  that  uploading  all  the  channels  will  take  an  unacceptable  length 
of  time.  It  may  be  useful  to  modify  the  control/simulation  program  to  prompt  the 
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user  for  which  channels  are  to  be  saved.  This  information  could  be  passed  to  the 
post-processor  through  the  .m  file. 
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5.  SYSTEM  MODELING  AND  CONTROLLER  DESIGN  FOR 

AFWAL  TRUSS 


The  truss  was  modeled  by  WL  and  the  model  was  given  to  OSU  in  the  standard 
F'EM  form,  with  no  reduction  of  model  order.  The  model  was  then  transformed  from 
nodal  coordinates  to  modal  coordinates  by  simultaneously  diagonalizing  the  mass  and 
stiffness  matrices.  The  model  order  was  then  reduced  using  balanced  model  reduction 
techniques  tis  described  in  Sections  2.2.1  and  2.2.2.  Controllers  were  then  designed 
in  continuous  time  for  the  reduced  order  models  then  discretized  for  implementation 
on  the  control  computer.  Designing  the  controllers  in  continuous  time  was  desirable 
since  the  sampling  rate  was  not  known  at  the  time  the  controllers  were  designed. 
Also,  many  of  the  design  tools  which  were  used  were  for  continuous  time  systems.  It 
appears  that  the  designs  did  not  suffer  much  from  this,  since  the  sampling  rates  were 
sufficiently  high. 

5.1  System  Modeling 

The  FEM  model  of  the  truss  was  constructed  using  MSC/pal  2  software.  The  model 
finally  used  was  an  equivalent  beam  representation  of  the  truss.  The  stiffness  prop¬ 
erties  of  the  beam  were  derived  from  a  single  bay  of  the  truss.  Using  a  single  bay 
limited  the  number  of  degrees  of  freedom  yet  allowed  an  accurate  approximation  of 
the  behavior  of  the  truss  for  the  lower  frequency  modes  of  interest.  The  beam  model 
divides  the  truss  into  16  elements  and  was  not  difficult  to  match  with  open  loop 
testing  of  the  truss.  The  full  model  has  33  nodes  and  136  degrees  of  freedom. 

The  eight  control  actuators  and  single  disturbance  actuator  are  modeled  as  simple 
springs,  masses  and  dampers,  as  shown  in  Figure  5.1.  They  are  constrained  to  move 
in  a  single  direction  with  respect  to  the  reference  frame.  In  actuality,  they  should 
be  constrained  to  move  in  one  direction  with  respect  to  the  truss;  however,  for  small 
deviations  in  the  truss,  it  can  be  assumed  that  this  is  equivalent. 

5.1.1  Modeling  of  Control  Forces  and  Sensor  Signals 

Control  forces  were  added  to  the  FEM  model  considering  the  action  of  the  proof  mass 
actuator.  The  force  produced  is  not  with  respect  to  a  stationary  reference  frame,  as 
would  be  the  case  with  a  thruster  or  an  actuator  attached  to  the  wall.  A  proof  mass 
actuator  produces  an  equal  and  opposite  force  on  the  proof  mass  and  the  body  of  the 
actuator,  as  indicated  in  Figure  5.1.  Recall  the  FEM  equation  (2.20)  is 

Mz  Dz  +  Kz  =  Q  (5.1) 
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Truss 


Mass 


Figure  5.1:  Model  of  Proof  Mass  Actuator 

where 

z  —  nodal  states  (displacements  and  velocities  in  each  direction,  at  each  node) 
M  —  mass  matrix  (symmetric) 

D  —  damping  matrix  (symmetric) 

K  —  stiffness  matrix  (symmetric) 

Q  —  generalized  external  forces  applied  to  each  node,  in  each  direction 

Define  the  input  to  the  actuators  to  be  u  and  the  output  of  the  velocity  sensors  to  be 
y,  where  the  sensors  are  numbered  as  shown  in  Figure  5.2.  The  u  and  y  are  directly 
related  to  the  nodal  states, 

y  =  C,z  (5.2) 

and 

Q  =  B,u  (5.3) 

Cg  has  as  many  rows  as  there  are  outputs,  and  Bg  as  many  columns  as  inputs. 
For  point  sensors  that  senses  velocity  along  a  single  coordinate  axis,  each  row  of 
Cg  contains  a  1  at  the  nodal  state  where  the  sensor  is  located,  and  the  rest  of  the 
elements  in  the  row  are  zero.  For  momentum  exchange  actuators  mounted  along 
a  single  coordinate  axis,  each  column  contains  a  +1  at  the  nodal  state  where  the 
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Figure  5.2;  Sensor  and  Actuator  Numbering  for  Controller  Designs 

actuator  is  attached  the  truss,  and  a  -1  at  the  nodal  state  of  the  proof  mass.  If 
relative  velocity  sensors  are  used  (such  as  the  LVT  sensors),  it  is  clear  that  C,  =  Bj. 


5.1.2  Transformation  from  Nodal  Coordinates  to  Modal  Coordinates 


The  MSC/pal  2  software  produces  a  FEM  model  with  272  nodal  states.  This  model 
is  transformed  from  nodal  coordinates  to  modal  coordinates  as  discussed  in  Sec¬ 
tion  2.2.1. 

The  FEM  model  produced  by  MSC/pal  2  software  contains  10%  damping  in  each 
proof  mass  actuator  mode  due  to  LVT  feedback,  but  includes  no  structural  damping. 
The  structural  damping  is  very  light  and  hard  to  model.  Dampings  of  3%  for  the  first 
bending  modes  and  1%  for  all  other  structural  modes  are  assumed.  This  damping  is 
added  directly  to  the  model  in  modal  coordinates  by  letting 


A  = 


0  / 

-n*  -^^D^-2c,n 


(5.4) 


where  is  the  diagonal  matrix  of  damping  ratios  and  0  is  the  diagonal  matrix  of 
squared  modal  frequencies. 


63 


Table  5.1:  Structural  Modes  and  Dampings  —  Open  Loop,  10%  Actuator  Damping 


Controlled 

1  Uncontrolled  I 

Mode 

1st  X 

1st  y 

1st  T 

2nd  X 

2nd  y 

3rd  X 

3rd  y 

2ndT 

4th  X 

4th  y 

3rd  T 

Frequency 

1.75 

1.75 

6.64 

8.46 

8.45 

19.42 

19.40 

20.22 

29.35 

29.31 

33.90 

%  Damping 

1.70 

1.70 

ES3I 

0.12 

0.12 

0.04 

0.04 

0.04 

0.02 

0.02 

0.02 

Table  5.2:  Structural  Modes  and  Dampings  —  Open  Loop,  50%  Actuator  Damping 


Controlled 

1  Uncontrolled  I 

Mode 

1st  X 

1st  y 

1st  T 

2nd  X 

2nd  y 

3rd  X 

liiiil!! 

2ndT 

4th  X 

4th  y 

3rdT 

Frequency 

1.67 

6.64 

8.44 

8.45 

19.39 

19.42 

20.22 

29.31 

29.35 

33.90 

6.32 

6.33 

1.08 

1.07 

1.07 

0.67 

0.68 

0.69 

0.62 

■il.-V 

0.61 

5.1.3  Model  Order  Reduction 

The  modal  transformation  described  in  the  last  section  results  in  a  system  equal 
in  order  to  the  FEM  model,  i.e.  a  model  containing  136  modes.  As  predicted  in 
Section  2.2.2,  this  model  has  too  many  modes  to  be  used  for  controller  design  and  the 
higher  order  modes  are  inaccurate.  Therefore,  the  model  was  reduced  as  discussed 
in  Section  2.2.2.  Two  models  were  created,  one  containing  14  states  and  another 
containing  42  states. 

The  controller  designs  for  the  truss  were  carried  out  using  a  14-state  model,  including 
5  structural  modes  and  2  actuator  modes.  A  common  problem  with  reduced  order 
controller  designs  is  the  destabilization  of  uncontrolled  modes.  To  test  for  this  con¬ 
dition,  the  42  state  model  was  used  <is  a  higher  order  “truth”  model  to  calculate 
closed-loop  damping  and  to  simulate  the  response  of  the  truss.  This  model  included 
1 1  structu.'al  modes  and  all  10  actuator  modes  (8  control  actuator  modes  and  2  modes 
to  model  the  disturbance  actuator  which  is  oriented  at  a  45°  angle  to  the  2  bending 
directions).  The  frequencies  and  natural  dampings  of  the  structural  modes  included 
in  the  truth  model  are  shown  in  Table  5.1.  The  actuator  modes  are  all  approximately 
1  Hz,  with  10%  damping. 

Table  5.2  shows  the  frequencies  and  natural  dampings  of  the  structural  modes  in¬ 
cluded  in  the  truth  model  for  the  czise  of  50%  actuator  damping. 


5.1.4  Effects  of  Discretization  and  Sampling  Delay 

In  the  simulation  of  the  truss  response,  the  plant  was  discretized  and  the  discrete¬ 
time  response  calculated.  The  discretized  plant  model  and  discrete-time  simulations 
|)roducefl  very  similar  results  to  those  of  the  continuous-time  model  and  continuous¬ 
time  simulations.  Thus,  discretization  did  not  present  any  problems. 

However,  a  more  serious  problem  was  encounted  with  delay  in  the  feedback  loop. 
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The  Optima  3,  due  to  the  way  in  which  the  A/D  and  D/A  cards  function,  has  a 
two  sampling  period  delay  between  output  and  input.  This  delay  was  not  taken  into 
account  when  the  controllers  were  designed,  although  it  is  modeled  in  the  simulation 
routine.  Several  of  the  controllers  have  been  destabilized,  or  nearly  destabilized  by 
this  delay. 

5.2  Controller  Design 

5.2.1  Performance  Requirement — Tip  Pointing  Error 

Some  mecisure  of  merit  was  needed  to  judge  the  performance  of  the  controllers,  as  well 
tis  actually  design  them.  For  this  criterion,  we  chose  to  minimize  the  displacement  of 
a  point  at  the  tip  of  the  truss  (24  inches  from  the  center).  The  light  source  is  located 
at  this  point,  so  the  performance  may  be  directly  measured  during  experiments  with 
the  optical  sensor  as  a  Line-of-Sight  (LOS)  error.  The  distance  of  this  point  from  the 
enter  of  the  truss  was  chosen  to  include  the  torsional  modes  of  the  truss. 

The  error  term  p  can  be  written  in  terms  of  displacements  of  the  top  four  sensors 
(yi,!/2,y4  and  r/s)- 

=  [^(yi +y2)  +  y  sin45®(yi -t/2)]  +[^(y4  +  ys)  +  y  sin45°(y4  -  ys)  (5.5) 
This  can  be  rewritten  in  matrix  form  — 


=  ulQ'Vd 

where  yd  is  the  vector  of  displacements  at  the  sensors. 
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(5.7) 


/:/  =  2.932  /i'  =  0.507  ^  = -1-219 


(5.8) 
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Now,  let 


y  =  CdX,  (5.9) 

where  Cj  is  the  output  matrix  yielding  displacements  at  the  sensor  locations.  This 
implies 

=  x^CjQ'CdX  =  x^Qx  (5.10) 

This  LOS  error  term  was  used  in  the  design  of  optimal  controllers. 


5.2.2  Controller  Design  Techniques 

This  project  presented  the  opportunity  to  test  a  wide  variety  of  control  techniques. 
Most  of  these  were  based  on  the  Linear  Quadratic  Regulator.  The  cost  function 
minimized  in  these  techniques  was 

J  =  +  u^Ru]  dt  (5-11) 

or  some  variation  thereof.  Most  of  the  controllers  were  designed  using  the  LOS 
requirement  described  in  the  previous  section  Jis  part  of  the  cost  criterion.  The 
inputs  were  usually  weighted  equally.  Thus,  we  typically  chose  R  =  a/,  where  a 
is  a  proportionality  constant  used  to  control  the  ratio  of  the  relative  weighting  of 
the  states  to  weighting  of  the  inputs  (which  is  equivalent  to  the  performance/control 
effort  trade-off). 

It  should  be  noted  that  the  limiting  factor  for  closed-loop  stability  was  typically  travel 
of  the  proof-masses  on  the  control  actuators.  In  the  controller  designs,  the  gains  grew 
larger  as  the  states  were  weighted  more  heavily  (or  equivalently,  the  inputs  weighted 
less  heavily)  to  increase  performance.  The  gains  grew  larger  which  required  more  force 
and,  hence,  further  travel  of  the  actuators.  Eventually,  a  point  was  reached  where  the 
high  gains  demanded  the  proof-masses  of  the  control  actuators  to  exceed  the  maxi¬ 
mum  pos.sible  travel.  Due  to  the  travel  limitations  of  the  proof-mass,  nonlinearities 
were  introduced  which  resulted  in  instability. 

Because  the  amount  of  travel  required  to  create  a  given  force  decreases  as  the  fre¬ 
quency  increases,  more  control  authority  is  actually  available  for  modes  above  the  first 
bending  mode.  Rather  that  using  the  LOS  requirement  in  the  cost  function,  these 
higher-frequency  modes  could  be  weighted  more  heavily-thus  increasing  the  perfor¬ 
mance  in  the  controlled  higher-frequency  modes.  Although  a  few  initial  controllers 
were  designed  using  equally  weighted  states,  we  munly  used  the  LOS  requirement  in 
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the  cost  function  and  did  no  ad  hoc  weighting  in  our  designs.  However,  it  should  be 
noted  that  our  results  could  perhaps  be  improved  with  such  ad  hoc  weighting. 

All  of  the  following  design  techniques  have  been  discussed  thoroughly  in  Section  2.3 
and  in  other  sources  [96];  therefore,  only  a  short  description  will  be  given  here. 


Linear  Quadratic  Gaussian  (LQG) 

This  approach  is  the  standard  full  state  feedback,  using  a  Kalman  filter  to  estimate  the 
modal  states  [50].  The  noise  in  each  input  and  output  is  assumed  to  be  uncorrelated 
and  equal,  so  W  and  V  are  proportional  to  the  identity  matrix.  The  relative  levels 
of  these  noises  are  unknown,  so  the  relative  weight  is  adjusted  until  the  response  of 
the  observer  is  satisfactory.  This  is  achieved  with  iW  =  V. 

The  Riccati  equations  were  solved  using  MATLAB  software  [94].  If  the  weighting 
matrix  Q  was  positive  semi-definite,  it  was  usually  necessary  to  add  a  small  number 
to  the  diagonal  for  numerical  stability. 

Linear  Quadratic  Gaussian  with  Loop  Transfer  Recovery  (LTR) 

This  approach  combines  the  frequency  domain  and  state  space  techniques  for  a  min¬ 
imal  phase  system  and  is  very  similar  to  the  LQG  design  discussed  above. 


Centralized  Optimal  Output  Feedback  (COFB) 

This  design  simply  minimizes  the  cost  function  given  in  the  decentralized  case,  except 
there  is  no  partitioning  into  subsystems.  Again,  DOLORES  was  used  except  in  this 
case  it  was  allowed  to  vary  all  the  terms  in  the  gain  matrix.  The  feedback  is  of  the 
form  u  =  Ky,  where  K  is  not  block  diagonal. 


Decentralized  Optimal  Output  Feedback  (DOFB) 

This  decentralized  design  is  based  on  the  following  four  groups  of  actuator /sensor 
(input/output)  pairs:  {1,2,3},  {4,5,6},  {7}  and  {8}.  The  feedback  is  u  =  Ky  where 
A'  is  8  X  8,  u  is  the  control  vector  and  y  is  the  vector  of  sensor  velocities.  K  is  block 
diagonal,  where  the  blocks  are  3  x  3,  3  x  3,  1  x  1  and  1x1.  This  design,  like  all  the 
other  output  feedback  designs,  was  obtained  with  software  called  DOLORES  which 
was  developed  at  OSU  to  solve  the  optimal  output  feedback  problem  [58,  97]. 
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Decentralized  Frequency  Shaping  (DFS) 

This  design  uses  the  same  decentralized  groups  as  the  regular  decentralized  case. 
The  cost  function,  however,  consists  of  the  plant  outputs  rather  than  the  states  [59, 
96,  98].  The  weighting  of  the  control  outputs  was  related  to  frequency,  in  order  to 
discourage  the  controller  from  acting  at  high  frequencies.  The  DOLORES  software 
also  incorporates  a  feature  to  augment  states  to  the  system  and  design  these  frequency 
shaping  controllers.  For  first  order  filters,  the  resulting  controller  is  an  8th  order 
dynamic  controller  of  the  form  u  =  K\y  +  K2Z  where  z  =  Fz  +  Gu.  K\  and  are 
block  diagonal,  partitioned  as  in  the  decentralized  case,  and  F  and  G  are  diagonal. 

Overlapping  Decomposition  (OD) 

This  method  divides  the  truss  into  six  disjoint  subsystems,  with  the  actuator/sensor 
(input/output)  pairs  grouped  as  {1,2},  {3},  {4,5},  {6},  {7}  and  {8}.  Each  subsystem 
was  minimized  separately  using  DOLORES.  The  feedback  isu  =  Ky  where  A"  is  8  x  8, 
u  is  the  control  vector  and  y  is  the  vector  of  sensor  velocities.  K  is  block  diagonal, 
where  the  blocks  are  2x2,  1x1, 2x2,  1x1, 1x1  and  1x1. 

Controlled  Component  Synthesis  (CCS) 

For  CCS  controller  design,  the  full  136  degree  of  freedom  FEM  of  the  truss  was 
decomposed  into  two  component  models.  The  lower  10  bays  comprised  the  lower 
component  which  had  62  degrees  of  freedom.  The  upper  6  bays  comprised  the  upper 
component  which  had  80  degrees  of  freedom.  These  component  models  were  formed 
from  the  composite  system  model  using  the  procedure  discussed  in  Section  2.3.6  since 
the  composite  model  was  available. 

Modal  decomposition  wais  performed  on  the  FEM  of  each  component  by  solving  the 
appropriate  eigenvalue  problems  for  the  component  modes  and  modal  frequencies. 
For  numerical  reasons,  the  lower  and  upper  component  models  were  reduced  to  20 
and  44  degrees  of  freedom,  respectively,  using  frequency  truncation.  At  this  stage,  all 
component  modes  below  100  Hz  were  retained.  A  state  space  model  w^ls  then  formed 
using  the  modal  states.  Taking  the  sensor  velocity  measurements  as  the  outputs  and 
the  actuator  signals  as  the  inputs,  the  component  models  were  further  reduced  based 
on  observability  and  controllability  of  the  modes.  The  resulting  models  were  used  for 
controller  design. 

Static  optimal  output  feedback  gains  were  found  for  each  component  so  as  to  min¬ 
imize  motion  at  the  internal  boundary  components  (the  one-half  and  three-quarter 
stations).  The  motion  at  the  tip  was  also  weighted  to  meet  additional  performance 


68 


objectives.  A  quadratic  performance  criterion  as  in  Equation  (2.110)  was  used  for 
the  component  controller  designs.  The  resulting  static  output  feedback  control  law 
for  each  component  was  of  the  form 

u  =  Kcc*y-  (5.12) 

Optimal  Projection  (OP) 

This  approach  uses  stochastic  modeling  for  the  uncertainty  of  the  system  order  to 
improve  the  robustness  with  respect  to  the  unmodeled  higher  frequency  modes.  The 
controller  is  of  the  form  u  —  Kx  where  x  =  A^x  +  Fy.  The  order  of  the  controller  is 
predetermined — in  this  case  ten. 

Decoupled  Optimal  Projection  (DOP) 

This  method  is  a  decentralized  version  of  optimal  projection.  The  controller  is  essen¬ 
tially  the  same,  except  three  third  order  controllers  run  independently,  one  for  each 
bending  direction  and  a  third  for  the  torsion. 


Hoo  Gap  Metric  Design 

The  goal  in  this  method  is  to  shape  (/  —  PC)~^P  for  disturbance  attenuation  and 
C{I  —  PC')"*  for  shaping  the  control  action  required.*  In  this  case,  only  the  top  four 
actuators  are  used  for  control.  The  multivariable  design  is  a  centralized  design  with 
the  four  sensor /actuator  pairs  as  the  inputs  and  outputs  of  the  system. 

A  balanced  truncation  of  the  42-state  “truth-model”  yields  a  suitable  design  model. 
The  balancing  transformation  and  truncation  are  performed  with  the  algorithms  dis¬ 
cussed  previously.  The  resulting  design  model  from  heis  14  states  corresponding  to 
the  first  five  structural  modes  and  two  actuator  modes.  The  gap  calculated  between 
the  modal  truncated  model  and  the  balanced  truncated  model  was  found  to  be  0.02. 
The  difference  between  the  two  models  is  primarily  due  to  the  actuator  modes.  Since 
actuator  modes  need  to  be  considered  for  stability  of  the  system,  one  might  ques¬ 
tion  whether  or  not  two  modes  are  enough.  To  answer  this  question,  both  reduced 
models  were  compared  to  a  30  state  model  which  contained  the  same  five  structural 
modes  and  all  10  actuator  modes.  The  gap  between  either  reduced  model  and  the  30 
state  model  was  found  to  be  around  0.09.  This  value  was  judged  to  be  acceptable; 
therefore,  a  fourteen  state  model  seemed  appropriate  for  design. 

'This  portion  of  work  was  performed  by  Scott  Buddie. 
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The  weighting  function  Wb  =  kW  is  applied  to  each  channel.  For  the  constant  weight, 
W  is  set  equal  to  I  and  k  is  varied  to  see  what  happens  to  the  resulting  closed  loop 
system.  Initially,  k  was  set  to  /  to  get  an  initial  indication  of  the  closed  loop  responses. 
By  doing  this,  the  modes  were  shifted  by  what  appeared  to  be  an  excessive  amount. 
After  some  additional  trials,  an  acceptable  range  of  k  was  found.  The  controllers 
were  calculated  using  the  four-block  solution  of  Englehart  and  Smith  [89]  for  reasons 
stated  later.  Therefore,  bopt  was  found  from  the  four-block  solution. 

In  case  the  14-state  controller  could  not  be  implemented  at  the  desired  sampling 
rate,  controller  reduction  schemes  were  employed  to  generate  12-state  controllers. 
An  upper  bound  for  the  weighted  controller  gap  between  the  full  and  reduced  order 
compensators  can  be  found  by  taking  twice  the  sum  of  neglected  singular  values  of  the 
product  of  the  controllability  and  observability  matrices.  This  controller  gap  comes 
out  to  be  about  0.16  for  A:  =  0.5. 

Since  the  optimal  controller  algorithm  was  causing  the  truth  model  to  become  un¬ 
stable,  we  looked  at  introducing  a  first  order  weighting  function  to  reduce  control 
action  in  the  higher  frequency  modes.  In  this  multivariable  case,  we  will  introduce 
the  first  order  weight  into  each  of  the  channels;  thus,  creating  a  weighting  function 
with  4  states.  This  tells  us  that  a  full  order  controller  using  this  weighting  scheme  will 
consist  of  21  or  22  states  depending  on  the  algorithm  chosen.  In  this  case,  we  have 
decided  to  use  the  optimal  controller  algorithm  to  see  if  the  truth  model  could  now  be 
stabilized.  The  controller  is  reduced  to  15  states  for  simulation  and  implementation. 
Two  first  order  weights  are  chosen.  The  frequency  chosen  as  the  roll-off  frequency 
was  7.162  Hz  (45  rad/sec).  This  was  chosen  somewhat  arbitrarily  to  give  a  little  less 
weighting  to  the  torsional  mode  and  the  2nd  x  and  y  mode.  The  frequency  could  be 
shifted  depending  on  the  weighting  desired.  The  difference  in  the  first  order  weights 
is  due  to  a  constant  term. 

The  controllers  are  discretized  at  a  sampling  period  of  5  ms.  This  was  determined 
to  be  a  sufficient  amount  of  time  to  process  the  control  law  and  was  fast  enough  in 
terms  of  the  modes  of  the  system  that  were  looked  at.  In  the  implementation  of  these 
controllers,  there  is  a  two  sampling  period  delay  due  to  the  A/D  and  D/A  boards 
as  mentioned  in  Section  4.2.4.  This  delay  could  cause  our  closed  loop  system  to  go 
unstable.  The  delay  is  of  the  form  where  T  is  10  ms  and  u;  corresponds  to  a  modal 
frequency.  Therefore,  the  higher  the  modal  frequency,  the  greater  the  phase  shift  that 
occurs.  To  investigate  this,  the  exponential  is  modeled  by  the  function  The 

function  is  then  augmented  to  the  plant  and  the  eigenvalues  of  the  closed  loop  system 
are  compared  to  the  eigenvalues  of  the  plant  without  delay.  The  eigenvalues  in  all 
cases  were  hardly  affected  by  the  delay;  therefore,  controllers  were  designed  for  the 
plant — without  the  delay  being  considered.  If  the  delay  was  determined  to  cause  a 
problem,  then  one  could  augment  the  plant  with  a  delay  function  and  do  the  design 
based  on  the  augmented  plant. 
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5.2.3  Output  Feedback  and  Actuator  Damping 


Through  the  course  of  the  designing  controllers,  it  was  realized  that  there  is  a  strong 
relation  between  the  addition  of  structural  damping  and  actuator  damping.  As  more 
damping  is  added  to  the  structural  modes,  the  actuator  modes  become  more  lightly 
damped.  This  can  be  seen  in  the  root  locus  for  constant  feedbeick  of  a  single  sensor 
and  actuator  pair.  The  structural  and  actuator  modes  move  in  different  directions  as 
the  feedback  is  increased.  This  is  true  since  the  proof  mass  actuator  gives  the  proof 
mass  and  the  truss  equal  and  opposite  forces,  and  only  the  velocities  of  the  truss  are 
available  for  feedback. 

There  are  at  least  two  solutions.  One  solution  is  to  feed  back  the  velocities  of  the  proof 
masses  in  addition  to  the  velocities  of  the  truss,  or  perhaps  only  the  relative  velocities. 
This  could  be  done  using  the  LVT  sensors.  There  is  another  solution,  however,  which 
is  simply  to  increase  the  natural  damping  of  the  actuator  modes.  Much  more  damping 
can  be  added  to  the  structural  modes  before  the  actuator  modes  approach  the  ju  axis 
and,  thus,  instability. 
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6.  SIMULATION  AND  EXPERIMENTAL  RESULTS 


In  this  chapter,  results  of  the  experiments  on  the  truss  are  presented.  In  general,  we 
have  shown  that  the  experimental  behavior  of  the  truss  is  very  similar  to  the  simulated 
behavior  of  the  model,  both  in  open  and  closed  loop  testing.  Several  controllers  were 
tested  using  various  inputs.  Although  some  of  these  controllers  did  not  adequately 
stabilize  the  truss,  the  results  typically  coincide  with  predicted  behavior.  As  the  Hoo 
controller  designs  differ  considerably  from  the  LQR-based  designs,  the  Hoo  results 
will  be  presented  in  a  separate  section  from  the  other  results. 

6.1  Noise  Levels  of  Velocity  Signals 

The  level  of  noise  in  the  velocity  signals  presents  one  of  the  first  problems  encountered 
in  this  experiment.  As  mentioned  in  Chapter  3,  the  signals  from  each  accelerometer 
are  integrated  to  provide  velocity  measurements  using  analog  circuits  built  into  the 
current  drivers  for  each  actuator.  The  analog  circuit  has  a  cutoff  at  0.1  Hz,  and  zero 
gain  at  DC.  Therefore,  the  frequency  response  achieves  maximum  magnitude  at  0.1 
Hz  and  rolls  off  at  -20  dB/decade  as  the  frequency  increases.  The  noise  from  the 
accelerometers  is  assumed  to  be  white,  however,  the  integrators  tend  to  amplify  the 
low  frequency  noise  to  the  point  where  it  could  be  a  problem. 

Table  6.1  shows  typical  noise  levels.  These  data  were  gathered  through  the  Optima  3 
while  the  truss  Weis  at  rest  and  represent  the  A/D  noise,  noise  in  the  integrator,  and 
noise  in  the  acceleration  signal.  The  low  frequency  noise  does  not  seem  to  cause  signif¬ 
icant  problems  for  the  controllers.  Any  noise  which  passes  through  the  controllers  is 
attenuated  by  the  truss.  The  main  problem  which  is  caused  by  the  noise  is  difficulty 
in  metisuring  the  performance  of  the  truss,  especially  in  terms  of  the  steady-state 
velocity  RMS  values. 

Obviously,  reducing  the  noise  in  the  signals  would  improve  the  quality  of  the  data. 
Reducing  the  noise,  however,  would  have  required  redesigning  the  integrators,  which 
would  have  been  costly  and  undesirable.  Another  option  was  to  increase  the  dis¬ 
turbance  signal  to  the  truss.  This  was  tried  and  achieved  some  measure  of  success. 
Perhaps  the  best  long  term  solution  would  be  to  avoid  using  integrators  which  by 
nature  amplify  low  frequency  noise. 


Table  6.1:  Typical  RMS  Noise  Levels  of  Velocity  Signals  (inch/s) 


Sensor 

#1 

#2 

#3 

#4 

#5 

#6 

#7 

#8 

Velocity 

0.019 

0.014 

0.048 

0.018 

0.020 

0.030 

0.021 

0.035 
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Figure  6.1:  Bode  Plot  of  Shaker  with  Voltage  Drive 
6.2  Disturbance  Generation 

As  previously  discussed,  the  truss  can  be  excited  from  either  the  disturbance  actuator 
which  is  mounted  on  the  truss,  or  the  shaker  which  is  mounted  to  the  wall  and 
connected  to  the  tip  of  the  truss.  Each  of  these  disturbance  generators  presents 
advantages  and  disadvantages. 

The  disturbance  actuator  was  limited  in  both  the  travel  of  the  proof  mass  and  the 
maximum  force  generated.  The  natural  damping  of  the  disturbance  actuator  could  be 
increased  to  reduce  the  travel  of  the  proof  mass,  but  this  would  also  change  the  model 
of  the  truss  slightly.  Even  if  the  natural  damping  was  increased,  the  maximum  force 
which  the  actuator  could  achieve  would  still  pose  a  problem,  especially  with  random 
signals.  The  maximum  output  of  the  D/A  converter  was  ±10  V  which  converts  to 
±1  lb.  Thus,  the  RMS  output  of  a  random  signal  was  limited  to  less  that  0.4  lb.  The 
maximum  D/A  output  was  indeed  the  limiting  factor  in  determining  the  maximum 
force,  not  the  power  amp  or  the  actuator  coil.  This  problem  could  be  solved  by 
amplifying  the  analog  signal,  but  a  suitable  amplifier  could  not  be  found. 

The  second  disturbance  generating  alternative  is  the  shaker  mounted  to  the  wall 
near  the  tip  of  the  truss.  This  device  has  the  potential  to  impart  a  much  larger 
disturbance  force  to  the  truss,  although  the  travel  of  the  tip  of  the  truss  is  limited  to 
1  inch  peak-to-peak.  At  first,  a  voltage  source  was  used  to  drive  the  shaker.  This  cid 
not  perform  well  at  the  resonant  frequencies  of  the  truss.  Figure  6.1  shows  the  Bode 
plot  of  the  shaker,  from  the  voltage  command  to  the  output  force,  measured  by  the 
force  gauge  at  the  end  of  the  shaker  arm.  The  dips  in  the  transfer  function  at  the 
resonant  frequencies  are  caused  by  the  back  EMF  of  the  shaker.  Driving  the  shaker 
with  a  current  source  eliminated  this  problem.  The  current  driver  for  the  disturbance 
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Figure  6.2:  Bode  Plot  of  Shaker  with  Current  Drive 

actuator  was  used  to  drive  the  shaker,  and  this  appears  to  have  worked  moderately 
well,  as  shown  by  Figure  6.2.  The  response  at  low  frequencies,  however,  remains  a 
problem.  Generating  larger  forces  at  lower  frequencies  would  require  more  travel  of 
the  shaker  arm.  At  the  first  bending  mode,  we  can  see  that  the  output  of  the  shaker  is 
about  71%  of  the  output  at  higher  frequencies.  This  is  due  to  the  poor  low  frequency 
response  of  the  shaker.  In  theory,  this  response  should  be  flat.  This  shortcoming  is 
magnified  because  most  of  the  response  of  the  truss  is  due  to  the  first  bending  modes. 
This  problem  explains  why  some  of  the  data  do  not  match  the  simulations  well  as 
would  be  expected. 

Perhaps  the  low  frequency  problems  with  the  shaker  could  be  corrected  by  prefiltering 
the  disturbance  signal  with  a  filter  having  a  transfer  function  which  is  approximately 
equal  to  the  inverse  transfer  function  of  the  shaker  over  the  frequency  range  of  in¬ 
terest.  The  filter  could  be  constructed  by  curve  fitting  the  desired  transfer  function 
over  the  frequency  range  of  interest  and  using  several  MATLAB  tools  which  would 
make  this  design  ecisy.  For  example,  the  command  yulewalk  could  provide  a  good 
approximation.  Also,  a  model  of  the  shaker  could  be  determined  and  included  in  the 
simulations. 

Another  solution  to  this  problem  is  to  find  a  shaker  that  hcis  a  better  low  frequency 
response. 

As  previously  discussed  in  Chapter  4,  the  disturbance  signal  is  generated  in  the  con¬ 
trol/simulation  program.  Figure  6.3  shows  the  spectrum  of  this  signal  leaving  the 
D/.A  as  recorded  with  the  OnoSokki.  We  chose  to  bandlimit  the  signal  to  50  Hz 
which  allows  us  to  put  more  energy  into  the  frequencies  which  contribute  most  to 
actually  exciting  the  truss.  This  frequency  range  also  includes  several  uncontrolled 
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modes,  so  the  results  would  allow  the  controllers  to  be  tested  for  susceptability  to 
spillover. 


6.3  Bode  Plots 

Bode  plots  of  the  structure  can  be  generated  relatively  easily,  using  the  OnoSokki  dual 
channel  spectrum  frequency  analyzer.  The  OnoSokki  obtains  both  the  magnitude  and 
phase  of  the  transfer  function  by  taking  the  Fast  Fourier  Transform  (FFT)  of  the  input 
and  output,  and  then  dividing  each  frequency  point.  This  process  is  repeated  and 
averaged  to  obtain  a  smoother  plot.  The  lower  the  frequency  range,  the  longer  it 
takes  to  obtain  a  smooth  plot.  For  0-50  Hz,  a  good  plot  can  be  obtained  in  about  one 
minute.  For  0-10  Hz,  at  least  three  minutes  is  necessary.  This  presents  no  problem 
for  the  control  program,  so  long  as  the  data  are  not  stored  in  the  Optima  3  for  the 
entire  length  of  the  experiment. 

Bode  plots  are  also  easily  obtained  from  the  data  which  are  uploaded  from  the  Optima 
3  to  the  Sun  Workstation.  A  MATLAB  routine  is  used  which  implements  basically 
the  same  method  as  the  OnoSokki.  The  experiments  were  performed  for  more  than 
164  seconds  to  provide  16,384  data  points  per  channel  per  experiment.  These  time 
data  were  used  to  compute  4096-point  FFT’s  for  15  windows  with  50%  overlap.  A 
Hanning  window  was  used  when  computing  the  FFT’s  of  each  window.  The  frequency 
responses  were  calculated  by  taking  the  ratio  of  the  FFT  of  the  output  to  the  FFT 
of  the  input.  The  resulting  15  responses  were  then  averaged  together  to  create  a 
smoother  plot. 

Open  loop  plots  for  both  experimental  and  simulated  responses  are  shown  in  Fig¬ 
ure  6.4.  Agreement  between  the  experimental  and  simulated  results  will  be  discussed 
in  the  next  section.  Only  Bode  plots  for  the  first  sensor  output  with  respect  to  the 
disturbance  input  are  shown.  The  Bode  plots  of  other  sensor  outputs  are  very  sim¬ 
ilar.  Note  that  bending  modes  in  the  x  and  y  directions  are  at  virtually  the  same 
frequencies  due  to  the  symmetry  of  the  truss. 

Closed-loop  Bode  plots  of  the  truss  with  several  different  controllers  are  shown  in 
Figures  6.5  through  6.9.  Each  Bode  plot  shows  the  experimental  and  simulation 
magnitude  comparisons.  Figure  6.9  compares  the  closed  loop  frequency  responses 
with  the  open  loop  response  from  experimental  data.  The  differences  between  the 
closed  loop  responses  can  be  hard  to  see,  thus  only  a  few  are  included.  These  Bode 
plots  are  perhaps  not  the  best  way  to  compare  the  controllers  as  the  controllers  were 
designed  to  minimize  the  motion  of  the  point  light  source.  The  actual  cost  criteria 
which  were  minimized  for  the  various  controllers  are,  in  some  cases,  theoretically  and 
fundamentally  different  even  though  the  objective  of  each  controller  design  is  the 
same.  Regardless,  these  responses  do  provide  significant  insight. 
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Figure  6.4:  Comparison  of  Experimental  and  Simulation  Open  Loop  Frequency  Re¬ 
sponses,  Velocity  ^1,  Shaker  1-lb  RMS  Random 
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Figure  6.5:  Comparison  of  Experimental  and  Simulation  DOFB  Frequency  Responses, 
Velocity  ^^1,  Shaker  1-lb  RMS  Random 
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Figure  6.6:  Comparison  of  Experimental  and  Simulation  COFB  Frequency  Responses, 
Velocity  #1,  Shaker  1-lb  RMS  Random 


Figure  6.7:  Comparison  of  Experimental  and  Simulation  Overlapping  Decomposition 
Frequency  Responses,  Velocity  #1,  Shaker  1  lb  RMS  Random 
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Figure  6.8:  Comparison  of  Experimental  and  Simulation  Controlled  Component  Syn~ 
thesis  Frequency  Responses,  Velocity  #1,  Shaker  1-lb  RMS  Random 


Figure  6.9:  Comparison  of  Experimental  Frequency  Responses  for  Several  Controllers, 
Velocity  #1,  Shaker  1-lb  RMS  Random 
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Frequency  responses  of  the  light  source  position  from  the  disturbance  input  to  the 
LOS  sensors  have  also  been  generated  from  both  experimental  and  simulated  data. 
When  the  experimental  data  were  taken,  the  gain  of  the  amplifier  for  the  photo  array 
signals  from  the  optical  sensor  was  set  incorrectly.  Therefore,  the  amplitudes  of  the 
time  response  data  were  off  by  a  factor  of  10.  Knowing  the  source  of  this  error,  the 
data  were  post-processed  by  multiplying  the  time  signal  by  10  or,  equivalently,  by 
adding  20  dB  to  the  frequency  response  to  give  the  correct  results. 


0 


•  100 


.1201 - - - ^ - - - - - - - ^ - ■ - ■ - - - 1 

0  5  10  15  20  25  30  35  40  45  50 

Frequency  (Hz) 

Figure  6.10:  Comparison  of  Experimental  and  Simulation  Open- Loop  LOS  Frequency 
Responses,  LOS  x.  Shaker  1-lb  RMS  Random 
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Figure  6.11:  Comparison  of  Experimental  and  Simulation  DOFB  LOS  Frequency 
Responses,  LOS  x,  Shaker  l*lb  RMS  Random 
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Figure  6.12:  Comparison  of  Experimental  and  Simulation  COFB  LOS  Frequency 
Responses,  LOS  x.  Shaker  1-lb  RMS  Random 
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Figure  6.13:  Comparison  of  Experimental  and  Simulation  Overlapping  Decomposition 
LOS  Frequency  Responses,  LOS  x,  Shaker  1  lb  RMS  Random 
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Figure  6.14:  Comparison  of  Experimental  and  Simulation  Controlled  Component 
Synthesis  LOS  Frequency  Responses,  LOS  x,  Shaker  1-lb  RMS  Random 
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Figures  6.10  through  6.14  compare  the  corrected  experimental  frequency  responses  of 
the  LOS  position  in  the  x  direction  with  the  responses  from  simulation.  Note  that 
these  frequency  responses  match  very  closely  at  the  modal  frequencies.  Figure  6.15 
compares  the  closed  loop  frequency  responses  and  the  open  loop  response  from  sim¬ 
ulated  data.  The  simulated  data  were  used  because  the  responses  are  smoother  than 
the  experimental  responses  and  thus  show  the  added  damping  more  clearly.  From 
Figures  6.10  through  6.14,  one  can  see  that  the  added  damping  shown  here  is  es¬ 
sentially  equivalent  to  the  added  damping  from  the  experimental  data.  Again,  the  » 

differences  between  the  closed  loop  responses  can  be  hard  to  see,  thus  only  a  few  are 
included. 

One  problem  with  the  post-processing  technique  described  above  is  evident  in  these 
responses.  Although  the  post-processing  corrected  the  signal,  it  also  amplified  the 
noise  to  levels  the  noise  never  actually  achieved.  This  is  evident  in  the  experimental 
responses  where  we  see  a  noise  floor  of  about  -60  dB.  It  is  assumed  the  actual  noise 
present  in  the  frequency  response  is  due  to  the  noise  in  the  photo  array  or  resolu¬ 
tion  (quantization)  of  the  photo  array  and  not  due  to  the  amplifier  itself.  If  this  is 
indeed  the  case,  then  had  the  amplifier  been  set  correctly,  the  noise  floor  would  have 
been  about  -80  dB.  Thus,  the  frequency  responses  generated  with  the  post-processed 
(corrected)  data  do  not  accurately  depict  the  true  noise  floor. 


6.4  RMS  Velocities 

The  link  between  the  Optima  .3  and  the  Sun  Workstation  is  very  slow  and  the  amount 
of  data  that  can  be  uploaded  is  therefore  limited.  The  RMS  values  of  the  sensor  and 
actuator  signals  are  very  useful  for  evaluating  controller  performance  and  can  be 
uploaded  in  less  than  ten  seconds. 

Table  6.2  shows  the  RMS  sensor  velocities  and  actuator  power  meeisured  experimen¬ 
tally,  with  a  1-lb  RMS  random  disturbance  from  the  shaker.  The  laist  column  is  the 
total  RMS  of  all  the  actuator  signals.  These  data  were  taken  for  164  seconds  after 
the  truss  was  allowed  to  reach  steady  state  for  30  seconds.  The  same  random  number 
generator  seed  was  used  for  all  experimental  trials. 

Table  6.3  shows  the  simulated  sensor  velocities  and  actuator  power  for  the  same 
controllers.  The  simulation  results  are  similar  to  the  experimental  results. 


6.5  Modal  Damping 

One  of  the  simplest  yet  perhaps  most  meaningful  way  to  evaluate  controller  for  a 
I..S.S  is  to  compute  the  closed-loop  damping.  This  provides  a  good  indication  of  the 
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Table  6.2:  Experimental  RMS  Sensor  Velocities  and  Actuator  Power  for  Random 
Disturbance,  1-lb  RMS 


type 

lOI 

#1 

#2 

#3 

#4 

#5 

#6 

#7 

#8 

Power 

Open 

ai 

0.295 

0.274 

0.198 

0.261 

0.261 

0.176 

0.168 

0.178 

—  "  ! 

DOFB 

Bl 

0.124 

0.122 

0.092 

0.125 

0.141 

0.135 

0.064 

0.080 

0.659 

COFB 

Kll 

0.130 

0.125 

0.105 

0.130 

0.136 

0.134 

0.077 

0.100 

0.596 

OD 

Dl 

0.137 

0.136 

0.123 

0.124 

0.142 

0.123 

0.088 

0.138 

0.408 

CCS 

Bl 

0.144 

0.136 

0.101 

0.164 

0.149 

0.146 

0.076 

0.106 

0.503 

Table  6.3:  Simulated  RMS  Sensor  Velocities  and  Actuator  Power  for  Random  Dis¬ 
turbance,  1-lb  RMS 


type 

lOI 

#1 

#2 

#3 

#4 

#5 

#6 

#7 

#8 

Power 

Open 

□1 

0.259 

0.268 

0.146 

0.260 

0.264 

0.139 

0.153 

0.156 

-J 

DOFB 

IBI 

0.138 

0.126 

0.092 

0.126 

0.141 

0.092 

BTiCTl 

0.080 

COFB 

IBI 

0.139 

0.127 

0.093 

0.129 

0.144 

0.098 

0.075 

0.079 

0.729 

OD 

0.154 

0.144 

0.102 

0.143 

0.155 

0.100 

0.101 

0.103 

0.517 

CCS 

IBI 

0.154 

0.150 

0.099 

0.157 

0.163 

0.101 

0.090 

0.098 

0.543 

performance  level  a  controller  will  achieve.  The  damping  ratios  for  each  controller 
were  calculated  using  the  42  state  model,  including  the  first  1 1  modes  as  well  «is  the 
10  actuator  modes.  This  was  done  in  continuous- time,  with  no  delay  in  the  feedback 
loop. 

For  each  type  of  controller  listed  in  Section  5.2,  several  controllers  were  designed  at 
various  times  throughout  the  project  with  various  weighting  matrices.  Every  type 
of  controller  designed  has  been  implemented  in  the  software,  however,  only  some  of 
them  are  presented  here. 

Controllers  were  designed  using  the  pointing  requirement  derived  in  Section  5.2.  The 
relative  weighting  of  the  Q  and  R  matrices  was  adjusted  so  the  damping  in  the  first 
mode  was  about  17%.  This  represents  the  most  damping  that  could  be  added  to  the 
truss  using  output  feedback  with  10%  actuator  damping.  If  the  relative  weighting  on 
the  matrix  R  is  reduced  any  further,  the  damping  remains  about  the  same  and  insta¬ 
bility  is  approached  ^ls  discussed  in  Section  5.2.2.  The  actuator  modes  become  very 
lightly  damped  which  increases  the  cost  faster  than  adding  damping  to  the  structural 
modes  decreases  the  cost.  The  closed  loop  damping  ratios  for  these  controllers  are 
shown  in  Table  6.4.  The  actuator  damping  column  is  the  lowest  damping  of  all  the 
actuator  modes. 
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Table  6.4:  Simulation  Closed  Loop  Damping  Ratios  (in  %)  —  Minimized  For  Pointing 


Control  Type 

1st  X 

1st  y 

1st  T 

3n<l  X 

2nd  y 

3rd  X 

3rd  y 

2ndT 

Actuators 

OL 

1.70 

1.70 

0.42 

0.42 

0.42 

0.34 

0.34 

0.34 

9.10 

DOFB 

11.49 

11.65 

6.10 

2.34 

2.35 

0.91 

0.91 

2.18 

3.35 

COFB 

10.13 

10.38 

6.78 

2.62 

2.56 

0.65 

0.65 

2.40 

4.07 

OD 

10.12 

10.38 

6.78 

2.62 

2.56 

0.65 

0.65 

2.40 

4.07 

CCS 

10.12 

10.38 

6.78 

2.62 

2.56 

0.65 

0.65 

2.40 

4.07 

Table  6.5:  Simulation  Closed  Loop  Damping  Ratios  —  Minimized  For  Pointing  Re- 


quirement,  50%  Actuator  Damping 


Control  Type 

KM 

1st  X 

1st  y 

1st  T 

2nd  X 

2nd  y 

3rd  X 

3rd  y 

2ndT 

Actuators 

OL 

1 

6.32 

6.32 

1.07 

1.07 

0.67 

0.68 

0.69 

49.35 

DOFB 

MM 

16.16 

17.13 

19.73 

4.88 

5.08 

1.61 

1.68 

6.56 

50.17 

COFB 

WBM 

17.24 

18.59 

16.34 

6.20 

6.79 

1.36 

1.48 

5.39 

47.40 

OD 

4 

18.09 

18.13 

17.13 

2.72 

2.72 

1.23 

1.23 

5.76 

48.91 

CCS 

14.76 

8.12 

3.03 

3.58 

1.28 

1.49 

3.03 

47.89 

A  second  group  of  controllers  were  designed  which  were  intended  to  illustrate  how 
additional  damping  could  be  achieved  by  increasing  the  natural  damping  of  the  ac¬ 
tuator  modes.  Damping  can  be  added  to  the  actuators  by  increasing  the  gains  of  the 
LVT  feedback  loops,  which  are  adjusted  by  a  10  turn  calibrated  dial  of  the  front  panel 
of  the  actuator  power  drivers.  An  actuator  damping  of  50%  was  chosen.  The  open 
loop  structural  damping  with  50%  actuator  damping  changes,  since  there  is  coupling 
between  the  actuator  modes  and  the  structural  modes.  As  mentioned  above,  the 
(lamping  ratios  of  the  closed-loop  actuator  modes  decrease  eis  the  damping  ratios  of 
the  structural  modes  increase.  By  adding  more  damping  to  the  open-loop  actuator 
inodes,  more  damping  may  be  added  to  structural  modes  of  the  closed-loop  system. 
The  closed  loop  damping  ratios  for  these  controllers  are  shown  in  Table  6.5  and  may 
be  compared  with  the  open  loop  damping  ratios. 

Attempts  were  made  to  estimate  the  experimental  damping  ratios  from  the  transfer 
function  data  using  a  circle-fit  damping  estimation  technique.  Unfortunately,  the 
results  were  very  sporadic  and  varied  tremendously  for  various  data  “windows”  to 
which  circles  were  fitted.  For  example,  when  estimating  the  damping  of  the  first 
bending  mode  for  the  DOFB  controller,  two  separate  attempts  from  the  same  trans¬ 
fer  function  data  (different  windows)  yielded  estimates  of  12%  and  36%  damping. 
These  variations  between  data  windows  may  be  related  to  how  well  the  circles  fit  the 
data.  With  such  wide  variation  in  the  results,  these  experimental  damping  estimation 
attempts  were  abandoned. 
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Figure  6.16;  Typical  Sinusoidal  Disturbance  Signal,  0.2-lb  RMS  1.75  Hz 
6.6  Time  Responses 

Time  responses  are  well  suited  for  examining  the  transient  responses  of  the  truss. 
Initial  responses  were  taken  by  applying  the  disturbance  and  control  simultaneously. 
Decay  responses  were  obtained  by  disturbing  the  controlled  truss  for  30  seconds,  then 
removing  the  disturbance  while  the  control  remained  on. 

The  damping  of  each  mode  can  be  examined  by  applying  a  single  frequency  sinusoid 
disturbance  to  the  truss.  Figure  6.16  shows  the  1.75  Hz  signal  used  in  the  initial 
responses  (Figures  6.17  and  6.18)  and  the  decay  responses  (Figures  6.19  through  6.20). 
The  time  constants  of  the  exponential  rising  and  falling  envelopes  do  not  appear  to 
be  the  same.  The  reason  for  this  is  not  known.  It  may  involve  friction  in  the  actuator. 
Note  that  these  responses  were  generated  with  10%  actuator  damping  and  controllers 
designed  for  that  configuration.  Therefore,  the  attenuation  and  damping  ratios  do 
not  match  with  the  discussions  from  previous  sections  where  actuator  damping  ratios 
were  50%. 


6.7  Hoc  Gap  Metric  Design  Results 

The  //oc  Gap  Metric  Design  experiments  addressed  different  issues  than  the  LQR- 
based  techniques  ais  discussed  in  Sections  2.3  and  5.2.2.  Therefore,  we  present  these 
results  separately.’ 

Table  6.6  shows  the  damping  ratios  of  the  controlled  modes  for  the  controllers  designed 


’This  portion  of  work  was  performed  by  Scott  Buddie. 
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Figure  6.17:  Initial  Response  -  Open  Loop,  Velocity  0.2-lb  RMS  1.75-Hz  Source 
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Figure  6.18:  Initial  Response  -  DOFB  #2,  Velocity  #1,  0.2-Ib  RMS  1.75-Hz  Source  1 
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Figure  6.19;  Decay  Response  -  Open  Loop,  Velocity  0.2-lb  RMS  1.75-Hz  Source 
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Figure  6.20:  Decay  Response  -  DOFB  #2,  Velocity  #1,  0.2-lb  RMS  1.75-Hz  Source  1 
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Table  6.6:  Damping  Ratios  of  Structural  Modes,  bgpu  and  Weighted  Gaps  for  Hoo 
Designs _ _ _ 


Mode 

k 

%  Damping 

kopt 

Weighted 

Gap 

1st  X 

1st  y 

IstT 

2nd  X 

2nd  y 

Frequency(Hz) 

1.66 

1.66 

6.64 

8.45 

8.44 

Open  Loop 

— 

6.3 

6.3 

1.1 

1.1 

1.1 

— 

— 

Constant 

Weight 

0.333 

16.14 

16.14 

2.78 

1.78 

1.78 

0.7756 

0.7761 

0.4 

18.93 

18.93 

3.26 

1.89 

1.89 

0.745 

0.82 

0.5 

23.21 

23.21 

3.99 

2.39 

2.39 

0.709 

0.88 

First  Order 
Weight 

0.33(0. 

15.80 

15.80 

2.18 

1.42 

1.42 

0.779 

0.40 

O.SOfd.li+45) 

.+45 

22.68 

22.68 

3.11 

1.76 

1.76 

0.712 

0.55 

Second  Order 

Weight 

0.50(0.1. +53)’ 
.’+26.5.  +  2809 

22.04 

22.06 

6.70 

4.60 

4.60 

0.699 

0.370 

0.65(0. 1.+52)-' 
.’+26.  +  2704 

25.53 

25.55 

8.55 

6.04 

6.04 

0.655 

0.446 

4 


with  the  given  weighting  function.  Only  the  modes  which  the  controller  was  designed 
to  control  are  shown.  Table  6.6  also  shows  the  calculated  values  of  bopt  and  the 
weighted  gap  between  the  design  model  and  the  truth  model  of  the  plant.  From 
this  table,  we  can  see  that  for  the  case  of  constant  weighting,  the  damping  ratios  of 
the  modes  increase,  the  optimal  robustness  radius  decreases,  and  the  weighted  gap 
between  the  design  model  and  the  truth  model  increases  as  the  constant  value  k  is 
increased.  At  A:  =  0.333,  the  values  of  the  weighted  gap  and  bopt  are  just  about  equal 
to  each  other. 

In  addition,  closed  loop  stability  of  the  truth  model  is  checked  and  verified.  This 
differs  from  the  case  of  the  optimal  controller.  The  damping  values  of  the  controlled 
modes  are  similar  using  the  optimal  controller;  however,  the  truth  model  wets  desta¬ 
bilized.  This  is  why  the  four-block  solution  is  employed  for  the  constant  weighting 
function  2is  discussed  in  Section  5.2.2. 

By  using  the  first  order  weighting  scheme,  the  optimal  robustness  radius  increzises 
very  slightly  and  the  weighted  plant  gap  decreases  by  a  considerable  amount.  This 
agrees  with  intuition  because  we  are  de-emphasizing  control  at  the  higher  frequencies. 
Note  that  the  damping  ratios  have  decreased  from  the  constant  weighting  case  in  the 
structural  modes  as  expected.  In  addition,  closed  loop  stability  was  checked  with  the 
truth  model.  For  these  controllers,  the  truth  model  v/as  stable. 

In  the  previous  designs,  a  large  increase  in  damping  is  evident  in  the  first  x  and  y 
bending  modes  without  much  of  an  increase  in  the  damping  of  the  torsional  mode  or 
second  bending  modes.  We  would  like  to  increase  the  damping  in  the  torsional  mode 
by  increasing  the  magnitude  of  the  weighting  function  in  the  frequency  range  around 
the  torsional  mode.  By  introducing  the  second  order  weighting  functions  shown  in 
Table  6.6,  the  damping  of  the  first  torsional  mode  and  the  second  bending  modes 
increase  because  they  are  penalized  more  heavily.  This  second  order  weight  is  applied 
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to  each  I/O  pair;  thus,  introducing  a  weighting  function  of  eight  states.  Therefore,  the 
optimal  controller  algorithm  yields  a  dynamic  controller  with  29  states.  Controller 
reduction  yields  a  21-state  controller  for  simulation  and  implementation.  Table  6.6 
shows  the  resulting  damping  ratios  and  the  values  of  bopt  and  the  weighted  gaps.  It 
should  be  noted  that  the  damping  ratios  for  the  first  six  uncontrolled  modes  were 
also  calculated  for  each  controller  design.  The  damping  ratios  for  these  uncontrolled 
modes  were  essentially  the  same  in  the  presence  of  the  controllers  as  they  were  for 
the  open  loop  case. 

Bode  plots  for  constant  weighting  with  k  —  0.5  and  for  the  two  second  order  weighting 
cases  are  shown  in  Figures  6.21  to  6.23.  These  plots  show  the  frequency  response  of 
the  first  sensor  output  velocity  with  respect  to  the  disturbance  input  force.  The  Bode 
plots  of  other  sensor  outputs  are  very  similar. 


Figure  6.21:  Comparison  of  Experimental  Open  Loop  and  Closed  Loop  Frequency 
Responses  for  Constant  Weight  Hoo  Design,  k  =  0.5,  Velocity  ^^1,  Shaker  1-lb  RMS 
Random 
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iirequcacy  (Ib) 

Figure  6.22:  Comparison  of  Experimental  Open  Loop  and  Closed  Loop  Frequency 
Responses  for  Second-Order  Weight  Hoa  Design,  k  =  '5*^28^’  Velocity  ^1,  Shaker 

1-lb  RMS  Random 


Figure  6.2.‘L  Comparison  of  Experimental  Open  Loop  and  Closed  Loop  Frequency 
Responses  for  Second-Order  Weight  Hrx,  Design,  k  =  Velocity  #1,  Shaker 

I -lb  RMS  Random 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  has  been  an  overview  of  a  fairly  extensive  study  on  applying  different 
control  techniques  on  a  cantilevered  truss  structure  for  active  vibration  damping.  A 
number  of  of  standard  issues  in  the  general  area  of  active  control  of  large  flexible 
structures  have  been  considered.  These  include 

1.  Modeling 

2.  Model  Reduction 

3.  Actuator  Dynamics 

4.  Performance  Specification 

5.  Implementation  Constraints  ' 

(a)  Computational  Delays 

(b)  Actuator  Locations 

(c)  Sensor  Locations 

(d)  Bandwidth  and  Authority  Limits 

(e)  Unmodeled  Dynamics 

(f)  Unknown  Parameters 

6.  Selection  of  Control  Approach 

7.  Computation  of  Feedback  Gains 

Again,  as  stated  above,  these  are  the  issues  that  any  control  designer  faced  with  a 
large  flexible  structure  will  encounter. 

On  the  other  hand,  the  truss  structure  under  consideration  was  a  fairly  simple  problem 
for  control  due  to  the  following  facts: 

•  The  behavior  of  the  structure  w«is  quite  similar  to  that  of  a  beam. 

•  Bending  vibrations  along  the  both  axes  and  the  torsional  vibrations  were  fairly 
uncoupled. 

•  The  vibrational  frequencies  (along  each  axis)  were  not  close  enough  to  create 
major  problems  in  controller  design. 
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These  simplifying  cispects  notwithstanding,  we  believe  the  experimental  structure 
with  the  multiple  actuators  and  sensors  and  the  computer  system  provides  a  very 
successful  testbed  and  furthermore  has  provided  both  the  OSU  Group  and  the  WL 
researchers  a  unique  learning  experience  that  is  invaluable  for  work  on  other  large 
flexible  structure  problems. 

As  discussed  in  the  preceding  pages  a  number  of  different  control  approaches  were 
analyzed  both  in  simulation  and  in  implementation  on  the  truss.  These  were: 

1.  Linear  Quadratic  Gaussian  (LQG) 

2.  Linear  Quadratic  Gaussian  with  Loop  Transfer  Recovery  (LQG/LTR) 

.3.  Centralized  Optimal  Output  Feedback 

1.  Decentralized  Optimal  Output  Feedback 

5.  Decentralized  Frequency  Shaping 

6.  Overlapping  Decomposition  (OD) 

7.  Controlled  Component  Synthesis  (CCS) 

S.  .Maximum  Entropy/Optimal  Projection  (MEOP) 

!).  Decoupled  Optimal  Projection  (DOP) 

10.  /Ax,  Gap  Metric  Design 

Each  one  of  the  approaches  implemented  Wcis  successful  in  vibration  damping.  It 
has  to  be  understood  that  control  design  approaches  are  based  on  specific  criteria 
such  as  minimizing  energy,  or  specifying  closed  loop  pole  locations  or,  reducing  effects 
of  parameter  variations,  or  rejecting  high  frequency  effects,  or  adapting  behavior  to 
emulate  specified  dynamics.  Thus  one  cannot  and  should  not  compare  these  designs 
on  a  single  criterion,  and  we  have  not  attempted  to  do  so.  Indeed,  we  have  attempted 
to  also  consider  approaches  that  have  computational  and  conceptual  advantages  when 
they  are  to  be  used  on  more  complex  structures  than  the  truss  configuration  here. 
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A.  Sample  Run  of  Control/Simulation  Program 


Integrated  acceleroaeter  or  LVT  feedback  (a/1)  ?  a 

Control  of  12b  Tmse  Response 
Ohio  State  University 

Niniana  saapling  tiaes  are  given  in  parenthesis 
Active  controllers  anst  be  run  at  this  saapling  tiae 

Actuator  daaping  at  SO*/. 

1)  Open  Loop  Response 

Designs  Miniaized  for  Pointing  Requireaent  -  Cheap  Control 

2)  Decentralized  Optimal  Output  Feedback  (2bs) 

3)  Centralized  Optiaal  Output  Feedback  (2as) 

4)  Overlapping  Decomposition  (2bs) 

5)  Controlled  Component  Synthesis  (2as) 

6)  LVT  Decentralized  Optimal  Output  Feedback  (2bs) 

7)  LVT  Centralized  Optimal  Output  Feedback  (2bs) 

8)  LVT  Overlapping  Decomposition  (2as) 

Designs  Minimized  for  Pointing  Requirement  -  Medium  Damping 

10)  Decentralized  Optiaal  Output  Feedback  (2bs) 

11)  Centralized  Optiaal  Output  Feedback  (2as) 

12)  Overlapping  Decomposition  (2as) 

13)  LVT  Decentralized  Optiaal  Output  Feedback  (2as) 

14)  LVT  Centralized  Optimal  Output  Feedback  (2bs) 

Designs  Minimized  for  Pointing  Requirement  -  Light  Damping 

15)  Decentralized  Optimal  Output  Feedback  (2bs) 

16)  Centralized  Optiaal  Output  Feedback  (2as) 

17)  Overlapping  Decoaposition  (2bs) 

18)  Controlled  Component  Synthesis  (2as) 

Designs  Based  on  Weighted  Gap  Metric 

20)  Full  Order  with  Constant  Weights  (Sas) 

21)  Reduced  Order  with  Constant  Weights  (Sas) 

22)  First  Order  Weight  -  Heavy  Damping 

23)  Second  Order  Weight 

24)  Smith  k  Georgiou  Constant  Weight  (Sas) 

25)  Smith  k  Georgiou  Second  Order  Light  (Sas) 

26)  Smith  k  Georgiou  Second  Order  Medina  (Sas) 

27)  Smith  k  Georgiou  Second  Order  Heavy  (Sas) 

Actuator  damping  at  10% 

30)  Open  Loop  Response 

Choose  type  of  control  (1-30)  1 
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Control  truss  or  siaulats  (c/s)  ?  s 
Enter  sasipling  rats  in  ■illisseonds  2 

Enter  data  recording  period  in  nillisecoads  (mltiple  of  sai^ling  rate)  10 
Disturbance  Sources 

1)  Disturbance  actuator  nounted  on  truss 

2)  Shaker  nounted  on  vail 


Choose  type  of  disturbance  excitation  (1*2)  2 

Loading  Model  Files 
getting  . /Batrices/c_42 . dat 
getting  ./Batrices/cd_42.dat 
getting  ./Batricss/ctravel42.dat 
Loading  Controller  Files  and  Discretizing 
Loading  Model  vith  SOX  Actuator  damping 
getting  ./Batrices/a_42_aS0_d2.dat 
getting  . /Batrices/bsH_42_aS0_d2.dat 
Enter  length  of  tine  in  seconds  for  data  recording  10 
Hunber  of  recorded  data  points  is  1000 
Is  this  OK  ?  (y/n)  y 

Enter  delay  in  seconds  before  beginning  recording  0 
Enter  length  of  tine  in  seconds  to  control  truss  11 


Disturbance  Signals 

1)  Bandwidth  Limited  Gaussian  White  Boise 

2)  Sinusoidal  at  one  frequency 

3)  Sinusoidal  Sweep 

4)  Burst  of  Bandwidth  Limited  Gaussian  White  Boise 

5)  Burst  of  Sinusoidal  at  one  frequency 


Which  type  of  disturbance?  (1-5)  1 
Enter  RMS  value  of  disturbance  1 
Cutoff  is  50Hz 

. . .  computing  disturbance  signal 

...  getting  ./disturbance_filters/fir50.2.dat 

corrected  mean  is  -0.000117 

corrected  std  is  1.000000 

Enter  filename  vith  no  extension  for  data  (will  be  put  in  ./data/  )  tmp 
Simulation  Loop  Running 
5*/.  Complete 
. . .  107.  Complete 
.  .  .  207.  Complete 
.  . .  507i  Complete 
. . .  BOX  Complete 
. . .  uploading  .n  File 
. . .  opening  data  file 

RMS  pointing  error  is  0.019642 
RMS  Force  is  1.010897 
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RMS  of  total  actuator  signal  is  0.000000 

RMS  of  actuator  *1  is  0.000000 

RMS  of  actuator  #2  is  0.000000 

RMS  of  actuator  #3  is  0.000000 

RMS  of  actuator  #4  is  0.000000 

RMS  of  actuator  #5  is  0.000000 

RMS  of  actuator  #6  is  0.000000 

RMS  of  actuator  *7  is  0.000000 

RMS  of  actuator  #8  is  0.000000 

RMS  of  vslocity  at  sensor  #1  is  0.248130  or  about  0.041355  volts 

RMS  of  velocity  at  sensor  #2  is  0.254297  or  about  0.042383  volts 

RMS  of  velocity  at  sensor  #3  is  0.138079  or  about  0.023013  volts 

RMS  of  velocity  at  sensor  #4  is  0.248013  or  about  0.041002  volts 

RMS  of  velocity  at  sensor  *5  is  0.253500  or  about  0.042250  volts 

RMS  of  velocity  at  sensor  #6  is  0.130227  or  about  0.021705  volts 

RMS  of  velocity  at  sensor  *7  is  0.153854  or  about  0.025642  volts 

RMS  of  velocity  at  sensor  #8  is  0.157053  or  about  0.026176  volts 

Maximun  travel  of  actuator  *1  is  0.040370 

Maximum  travel  of  actuator  #2  is  0.043281 

Maximum  travel  of  actuator  #3  is  0.026063 

Maximum  travel  of  actuator  *4  is  0.043110 

Maximum  travel  of  actuator  #5  is  0.040240 

Maximum  travel  of  actuator  #6  is  0.025659 

Maximum  travel  of  actuator  *7  is  0.014906 

Maximum  travel  of  actuator  #8  is  0.014914 

Maximum  travel  of  actuator  *9  is  0.055187 

•m  file  is  in  ./data/ 

. . .  uploading  data 
vriting  0  buffers:  . 
writing  0  buffers:  . 
writing  0  buffers:  . 
writing  0  buffers:  . 
writing  0  buffers:  . 
writing  0  buffers:  . 
writing  0  buffers:  . 
writing  0  buffers:  . 
writing  0  buffers:  . 
writing  0  buffers:  . 
writing  0  buffers:  . 

.mat  file  is  in  ./data/ 

done 
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